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PREFACE 


"The  wholci  thing  is  a  lov  put-up  job  on  our  noble  credulity,"  said  Satn. 


NORNLfLN  LINDSAY,  The  Mae,ic  Pudding 


if  you  can't  dazzle  'em  with  brilliance. 
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The  clectroa-phonon  interaction  in  solid  argon  and  neon  was 

calculated  ^/ithin  tb.e  self-ccnsisrent  field  5Uf:m.;;nted  plane  wave  local 

exchange  (SCF  APK-Xa)  model.  The  Frb'hlich  HamiLltonian,  which  uses  the 

felectron-nhonon  interaction  in  the  form,  of  a  trarisition  probabiJ.ity 

matrix  elenent ,  was  assumed  to  be  appropriate  for  describing  correctioxis 

to  the  phouon  spectra.  Second  order  perturbation  theory  involving  the 

calculated  (.led ron-phonou  matrix  element  (EPME)  was  used  to  compute 

correctioTii.  to  obonoTis  unperturbed  by  the  elcc tion-phonon  interaction. 

The  unperturbed  phonons  were  calculated  within  the  self-consistent 

harmonic  (SClO  approvim.ation ,  assunJaig  a  superpo.^ition  of  Lcm;ard-Jone3 

pair  potential  r.  Tlie   pair  potential  paraiueter;.-  were  not  determ.ined  empiri- 

callv,  but  v/ere  fitted  ro  the  calculated  equilibrium  lattice  constant  or 

the  experim.enta2  lattice  constant  and  to  the  corresponding  calculated 

APV!-Xa  cohesive  energies.  Explicit  expressions  for  the  EPME  were  derived 

in  the  APW  formalism,  as  Xs-ere  selection  rules  bcrvreen  representations 


r 


involved  in  the  relevant  inter-band  transitions.  To  calculate  the  EPME 
nunerically,  computer  codes  vere  written  and  tested.  The  perturbed  pho- 
nons,  calculated  as  a  test  of  the  EPME,  vrere  found  to  be  in  poor  agreenent 
for  neon  and  in  good  agreement  for  argon,  a  distinction  which  can  be  traced 
directly  to  the  shortcomings  in  the  original  Xa  calculation  on  neon,  and 
not  to  the  method  in  eeneral. 


CHAPTER  I 


INTRODUCTION 


1-1.  Purpo se  of  Ca Iculnticn 

It  is  essential  to  study  the  interaction  of  electrons  with  the  lat- 
tice vibrations  or  phonons  in  a  crystal  if  the  basic  properties  and 
characteristics  of  solids  are  to  be  understood.  Such  effects  as  super- 
conductivity in  r.ietals,  phonon  dispersion  anor.alieSj  and  electronic  speci- 
fic heat  enhancement  may  be  calculated  once  one  has  detailed   knowledge. 
of  the  electron-phonon  (EP)  interaction. 

However,  the  airc  in  this  work  is  not  so  ambitious  as  a  calcuj.ation 
of  the  above  effects.  Rather,  the  present  effort  is  to  develop  a  formal 
computational  procedure  to  calculate  from  first  principles  the  Ef  inter- 
action in  a  physically  realistic  nocel  of  an  insulator.  Once  the  method 
has  been  justified  and  demonstrated  to  be  feasible,  the  application  to 
tlie  effects  caused  b}'  the  EP  inceracti.oa  ray  be  aone . 

A  seeniinj,ly  practical  and  jound  mcthad  o,^  calculating  the  E?   inter- 
action wa?  advanced  alinost  10  years  ago  by  Golibersuch.   His  idea  for 
calcu.latinc  the  elect ron-pboron  matrix  element  (EPIIi',)  v;as  to  exploit 
certain  features  cf  Slater's  augmented  plane  wave  (APW) "  method  for  the 
calculation  of  the  eneigy  bar. us  in  a  r.oiJd. 

In  his  article,  GcliLerr;uch  justified  a  certain  specific  foxm  of  the 
EPME,  then  explicitly  derived  the  computational  forrf.ulae  in  terms  of  som.e 
oi  I  ne  qi..;iui  ii_.n--r.  Qoailauj.o  ii.  cocvenu.iunaJ.  .-.i  ..  co4i.putei  coae„.   Jiov-e .  ej. , 


no  calculations  based  on  his  method  have  ever  appeared  in  the  lite.rature. 
I'he  reason  for  this  absence  of  calculations  is  probably  the  complicated 
nature  of  his  formulae;  in  particular,  the  occurrence  of  doubJ.e  sums  over 
angular  monienta  (two  "•  sums  and  two  m  sums)  which  include  some  rather  in- 
tricate spherical  harm.onic  transformation  matrices. 

An  exam.inaticn  of  Golibersuch's  paper  shows  that  his  metliod  does  not 
take  explicit  advantage  of  the  muffin-tin  symmetry  inherent  in  the  usual 
APW  method."^  Also  he  does  not  consider  an  insulator,  although  his  justi- 
fication of  the  EPME  is  general.  Therefore,  a  method  that  exploits 
the  muffin-tin  symmetry  has  been  developed  for  filled  valence  band  soJ.ids 
(not  necessarily  insulators),  and  is  the  subject  of  the  present  work.  It 
should  be  emphasized  however,  that  the  EPf-IE  used  in  this  work  is  based  on 
the  form  justified  by  Golibersuch. 

To  test  the  calculation  of  the  EPME,  the  effect  of  the  electro-  "inter- 
action on  the  phonon  spectra  of  solid  argon  and  neon  has  been  coirputed. 
Though  this  test  is  rather  severe  (as  will  be  explained  subsequently),  it 
demonstrates  both  the  power  and  the  limitations  inherent  in  our  procedure. 

3.-2  ■  As  sump  tlons  and  Approxin:ations 

Perhaps  the  iriost  crucial  assumptica  for  the  "verall  calculation  is 
that  of  allowing  tb.e  nucHei  and  electrons  to  interact  among  them.se]  ves 
and  each  ether   only  through  short-ranged  screened  potentials.  Such  screened 
potentials  are  inlicrrnt  in  the  se]  f-consist  i.it  field  (SCF)  APW  description 
of  the  electroas  and  the  Frclhlich  lami  Ito;  J  sn''  which  will  be  used  to  calcu- 
late energy  corrections  to  the  phoiions  aricing  from  the  EP  interaction.  A 
truly  basic  approach  v;ould  be  to  calculate  first  the  bare  nuclear  inter- 
actica  and  then  screen  this  interaction  witli  an  eleccron  distribution 


whose  response  is  described  by  a  dielectric  function.  An  approach  along 
these  lines  has  been  derived  by  Sinha^  for  the  APW  method,  but  does  not 
appear  to  be  computationally  feasible.  Although  the  formal  theory  of 
microscopic   lattice  dynamics  has  been  developed  by  Toya,   Pick,   and 
others®  through  the  dielectric  response  concept,  few,  if  any,  practical 
numerical  calculations  have  been  attempted.  This  work,  even  with  its 
approxiiuatioas,  represents  an  attempt  at  a  calculation  based  in  principle 
on  microscopic  lattice  dynamics. 

The  Born-Oppenheimer  adiabatic  theory  and  approximation   are  assumed 
to  be  applicable  in  calculating  the  electronic  wave  functions  for  a  static 
lattice.  The  foircal  equivalence  between  the  Fryhlich  theory  and  the  adia- 
batic theory  has  been  pointed  out  previously,    but  will  be  derived  expli- 
citly in  Cb.apter  II  in  the  present  context. 

In  order  to  carry  out  an  APV-'  calculation  of  the  now-traditional 
muffin-tin  type  (the  sort  originally  discussed  by  Slater  ) ,  a  one-electron 
potential  must  be  available.  Even  at  the  level  of  Hartree-Fock  theory, 
the  generation  of  such  a  potential  for  a  crystal  is  an  almost  prohibitively 
difficult  task,  A  great  deal  of  effort  has  gone  Into  developing  methods 
vrhich  circu-ivent  this  difficulty  at  the  cost  of  the  introduction  of 
certain  ad  hoc  physical  assui-ptions .  The  t.iost  widely  studied  of  these 
schemes  approxiinaces  exchanf';e  and  correlation  contributions  to  the  poten- 
tial along  the  lines  of  exchange  (designated  as  X)  in  the  electron  gas 
with  the  added  flexibility  of  one  parameter  (designated  as  o:)  :  the  model 
is  therefore  kr.ov.-n  as  the  "Xu"  model. 

In  this  study  of  lattice  dynanics ,  the  APU  r.uf  fin-tin  potential 
associated  with  a  particular  nucleus  is  assumed  to  move  rigidly  with  that 


nuc 


leus.  This  "rigid  ion"  approximation-^  allows  the  total  solid  potential 


to  be  a  superposition  of  all  muffiii-tin  potentials  centered  on  the  nuclear 
locations. 

For  the  nuclear  motion,  the  potentials  are  assumed  to  be  hartronic; 
i.e_.  the  effect  of  the  force  is  linear  in  the  displacement.  Even  the  self- 
consistent  harmonic  (SCH) '   approAimation  used  to  calculate  the  unperturbed 
phonons  has  as  its  basis  functions  harmonic  solutions,  though  the  method 
itself  contains  implicit  and  inherent  multi-phonon  processes  (to  infinite 
order).  ^"^  The  total  static  lattice  potential  is  assumed  to  be  describable 
in  terms  of  the  cohesive  energy  and  lattice  constant  of  the  solid.  This 
description  is  parameterized  through  a  set  of  Lennard-Jones  pair  poten- 
tial constants."^  The  fitted  constants  are  used  in  a  Lennard-Jcnes  pair 
potential  which  is  superposed  to  yield  the  total  lattice  potential. 

All  these  assumptions  and  approximations  will  be  reviewed  and  examiiied 
in  the  proper  context  in  Chapter  II  where  the  formal  theory  is  developed. 
In  light  of  the  conclusions  reached  in  Chapter  IV,  the  severity  or  appro- 
priateness of  the  approximations  e.'xd   assumptions  w-ill  be  discussed. 

1-3.  Out line  of  Calculation 

The  first  phase  of  the  calculation  consists  in  obtaining  the  static 
electronic  Bloch  wave  functions  at  a  given  l^rttice  constant  in  the  /J^W 
basis.  The  niimber  of  eigenvalues  and  Eloch  wave  fxractions  are  determ.ined 
by  the  type  and  number  of  atoms  in  th.e  solid.  Self-consistency  is  obtained 
throi;gi)  one  of  three  criteria:  coiiveigence  (or  m.ore  exactly,  iteration  to 
iteration  stability)  on  the  pressure,  one-electron  eigenvalues  (i^.e_.  band 
energies).,  or  potential.  Usually  when  one  of  these  criteria  is  achieved, 
the  others  are  also  fairly  well  satisfied, 

ITaau    the  Kloch  fuictions,  e.iOrgy  eigenvalvics,  and  self ■  consistent 


potential  are  obtained  to  desired  accuracy,  the  cohesive  energy'  of  the 
solid  is  calculated  at  the  assunied  lattice  constant.  The  cohesive  energy 
and  the  assunied  lattice  constant  then  can  be  used  to  parameterize  the 
electronic  contribution  to  the  bulk,  static  solid  in  terms  of  Lennard- 
Jones  pair  potential  parameters. ^ ^  This  parameterization  describes  the 
static  solid  by  a  lattice  sum  involving  all  distinguishable  pairs  in  the 
solid.  The  12-5  form  of  the  potential^"  in  used  for  cor.putational  conve- 
nience in  obtaining  the  phoncri  spectrum  of  the  solid. 

Next,  the  SCH  formalism  is  applied  to  the  given  total  lattice  poten- 
tial, based  on  the  Lennard-Jories   potential  superposition.  A  phonon  spec- 
trum is  calculated  to  a  specified  level  of  precision  by  iterative  solution 
of  the  SCH  dynamical  matrix  eigenvalue,  problem  (see  Appendix  I)  .  This  set 
of  phonons  constitutes  the  unperturbed  pb.onons  (i_.ft^.  phonons  unperturbed 
by  explicit  interaction  with  the  electrons). 

In  principle,  the  EPMr]  should  be  calculated  for  all  cransitions  be- 
tv/een  occupied  bands  and  empty  bar'ri,-.  Since,  hov.-ever,  the  transition  prob- 
ability is  dominated  by  a  factor  iirversely  nroi. ortional  to  the  energy 
difference  between  the  states  of  the  transition,  the  main  contributing 
transitions  are  takei  to  be  belween  the  higiiest  single  valence  band  states 
and  the  lowest  single  conductio:i  bai-.d  stat-'s.  In  the  rare  gas  insulators 
(e.£.  argon  and  neon),  this  means  only  transitions  from  the  highest  p-like 
valence  band  st.Ttes  to  the  lowest  s-like  ccnducticn  band  states.  The  fact 
that  these  bands  exbibit  simple  syiimietry  properties  throughout  the  zone 
is  a  great  convenieiice  for  ccirpr  ts:tion . 

The  I>';iE  is  calculated  throughout  the  entire  first  Erillouin  zone 
(BZ)  (see  Appendix  II).  7'lie  square  of  the  EP'.IE  is  used  to  find  a  second 
urdii-  pei  LUi-baLiD:i  ci.crgv  coi';:>.-ctic,a  to  the  pIirjnOi\s.  The  contrJ 'M.K:ions 


to  the  energy  shifts  for  a  pbotion  of  a  given  vector  and  polarization 
are  governed  by  certain  selection  rules  that  relate  this  phonon  to  the 
electronic  irreducible  representations  (which  label  the  electron  v'ave  vec- 
tor and  state) .  These  allowed  electronic  corrections  are  then  suiraned  over 
the  entire  first  BZ . 

This  is  a  basic  description  of  the  procedure  involved  in  the  EPME 
and  phonon  correction  calculation.  Chapter  II  gives  a  theoretical  and 
formal  justification  for  the  procedures  outlined  here.  Chapter  III 
includes  the  results,  coininentary,  and  comparison  with  known  quantities. 
Conclusions  regarding  the  results,  their  implications  for  the  assumptions, 
and  the  method  in  general  are  presented  in  Chapter  IV.  There  are  four 
Appendices  in  vrhich  are  explained  the  potential  and  SCH  approxim.atlons 
(Appendix  I),  specific  group  theory  and  selection  rules  (Appendix  II), 
the  actual  formulae  for  ti  i'  EPME  (Appendix  III) ,  and  a  description  of  the 
computer  codes  developed  and  used  (Appendix  IV) . 


CHAPTER  II 

FORMAL  THEORY  AIN'D  METHODS  OF  CALCULATION 

2-1.  Bora-Op penheimer  Tht^ory  for  Crystalline  Solids 

Crystalline  solids,  by  definition,  possess  a  definite  periodic  trans- 
latioaal  syn-jnefr;/.  Though  the  specific  crystal  structure  is  determined  by 
minimization  of  the  free  energy,  the  various  lattice  s^'Ttime tries  inherent 
in  a  particular  crystal  structure  frequently  may  be  exploited  to  reduce 
the  difficulty  in  carrying  through  whatever  formal  ansatz  has  been  c].ected. 
Thus,  the  translati.onal  periodicity  of  a  perfect  crystal  aj.lows  one  to 
reduce  manj'  problems  to  considerations  vrithin  a  single  cell  by  application 
of  Bloch's  theorem.  Group  theoretical  methods  also  may  be  used  to  take 
exp].icit  advantage  of  tlte  crystal  symmetry;  as,  for  example,  in  the  b^i  ock 
diagonal! zation  of  secular  matrices  ccnmonly  done  in  energy  band  and 
molecular  theory. 

The  eleaientary  form  of  Eloch's  theorem  for  a  particle  moving  in  a 
p■^ri.odic  potential  V(:^)  is 


y,  (r+R  )  -   exn[ik.K  ]':'  (r)  (2-1) 

k  ~-  ■  u  n  jc  - 

V(r4R  )  =-  V(r)  (2-la) 


vhere  \"  (r)  is  the  particle  e:i  genf  unction  with  v.'ave  vector  Ic  and  energy 
k 

E  ,  and  R   is  any  trauslaticn  that  leavris  the  solid  or  lattice  invariant, 
k      -n 

The  extension  to  m.ulti-particlc  systems  is  straightforward. 


The  symmetries  characterising  a  solid  are  valid  only  for  a  static 
crystal  structure.  In  reality,  the  solid  coupists  of  moving  nuclei  v/hose 
motion  is  centered  on  the  lattice,  sites.  P.owe-"er,  because  the  nuclei  are 
bound  in  the  crystal,  they  appear  on  the  average  to  be  at  the  equilibrium 
(i^.e_.  static)  lattice  sites.  It  is  then  inr.uitively  plausible  (because  of 
the  smallness  of  the  electron  mass  relative  to  the  nuclear  mass)  to  begin  a 
treatment  of  electrons  in  a  solid  by  assuming  that  the  nuclei  may  be  taken 
as  fixed  at  the  given  lattice  sites  for  a  specified  density.  After  a  solu- 
tion to  this  problem  is  achieved,  the  motion  of  the  nuclei  may  be  taken 
into  account  by  a  perturbative  calculation  of  the  electronic  response. 
Just  such  an  endeavor  is  the  purpose  of  this  research. 

The  formal  basis  for  this  kind  of  effort  was  first  enunciated  by  Eorn 
and  Oppeuheimer^  and  rhe  general  theory  has  become  known  by  their  names. 
Folloving  the  notation  of  references  10,  16,  and  17,  the  total  solid 
Hamiltonian  H  is  written  as 

H(R)  -   H  (R)  +  H  (r,R) 
—      n       e  ^  ^ 

=  0    P^/2M  +  V  (R)  )  +  (T  p^  +  V   (r,R)  +  V  (r)  )   (2-2) 

^.2  n  -       ;;   1    ne  e  ~ 

1  i 

where  Rydberg  atomic  units  v;ill  be  employer'  throughout  this  work  except 

as  othe:-t.'ise  noted.  In  (2-2),  M  is  the  nuclear  mass  v.-ith  the  sum  j  over 

nncl'^i,  tht->  sum  i  is  over  electrons,  V  (V)    is  the  nuclear-nuclear  inter- 

n  -~ 

action  potential,  V   (r,r.)  is  che  nxicleiT-eL^ciron  interaction  potential, 
ne 

and  V  (r)  is  the  electron-electron  interacta-^n  potential.  Also  R  is  a 
e  ~" 

supervector  of  all  nuclear  position  vectors  and  £  plays  a  similar  role 
for  electrons.  The  grouping  of  the  terr.iS  in  (2-2)  is  suggestive  of  the 
role  v.'hich  the  nuclear  coordinates  play  as  parar^ieters  in  the  adiabatic 
or  instant  anc-ons  lesnonse  theory  of  el'ictronl;:  structure.  That  is,  the 


last  threo  terms  in  (2-2)    give  the  Schrt'dingei  eqn£\tion  for  the  electronic 
state  as 

H  (R)^',  (r,R)  =  6   p^  +  V   (r,R)  +  V  (r)  )T,  (r,R) 
1  — 

=  E^(R)H'^(r,R)  (2-3) 

th 
with  E,_(R)  the  total  electronic  energy  of  the  k  *  state  dependent  upon 

the  nuclear  coordinates  R  as  paraiaetcrs. 

The  electronic  eigenf unctions  of  (2-3)  forra  in  general  a  complete 

basis  set  and  therefore  may  be  used  to  expand  the  total  system  (nuclear 

and  electronic)  wave  function 

r.(r,R)  =  y  'Hk,R)'t',  (r,R)  (2-4) 

k        -- 

where  the  expanse  en  coefficients  depend  on  the  k  '  electronic  state  and 

the  miclear  coordinates  R.  Operating  on  (2-4)  with  the  total  Hamiltonian 

H  (2-2)  gives  three  tei'ms 

KH(r,R)  =  );(  \(r,R)[y  p^/2M  +  V^(R)  +  Ej^(R)  ]<;  (k,R) 
'  '     k   -      j   -^       ^       - 

-  1/M  V^\(r,lO-Vj,4>(h,R)  -  1/2M  *  (k,R)V^'I'j^(r  ,R)  )   (2-5) 

By  multiplying  the  last  two  terms  in  (2-5)  by  'ij;,(il,R)  and  integrating 
over  the  electronic  coordinates,  these  terms  may  be  shown     '  to  be 
related  to  electronic  scattering  fror.i  the  state  k  to  k'.    Because  this 
scattering  is  not  in  accordance  v.'ith  the  adiabatic  theory  (which  requires 
the  electronic  state  be  unchanged)  ar.d  because  the  terris  are  small  due  to 
the  factor  of  the  nuclear  mass,  the  last  tvjo  terms  may  be  eliminated  from 
coasideration . 

Examination  of  the  rest  of  (2-5)  (after  integration  over  the  elec- 
trcni,-  ccordinalcs  and  elimiaatica  cf  the  n'T  -  ad^ -hatic  tf'r'^is)  she'-s  that 
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ve  have,  achieved  a  quite  plausible  separation  between  the  nuclear  and 
the  electronic  motions.  This  separation  is  emphasized  by  requiring  that 
the  <l-'(k,R)  be  eigsuf unctions  of  the  nuclear  motion  determined  from  (2-5) 

(  E^(_R)  +  E^.(_PO  )  ^   Ak,K>    =  E   (R)$   (k,R)     (2-6) 
where,  ia  ^.he  case  of  a  lattice  dynamical  treatment, 

\}S^^   -'    \X   -^  1/2)  .^  (2-7) 

E  ,  is  the  total  energy  of  the  n    state  for  the  £_  '  phonon  of  frequency 
Uj   and  A  polarir,ai  ion  . 

The  separation  of  the  electronic  and  nuclear  motion  has  been  accom- 
plished by  use  of  the  adiabatic  hyxioc.liesis.  The  nuclear  motion  governed 
by  (2-6)  raay  now  be  described  as  a  ph.oncn  field  either  of  the  ordinary 
harL^:cnic  typ<;,  or  the  intrinslcaj.ly   quanrun  mechanical  type  given  by 
SCH  theory  (see  Appendix  I) .  Roth  approaches  diagonalize  a  force  constant 
or  dynamical  matrix  to  obtain  the  phonon  eigenvalues.  The  only  difference 
is  in  the  detei-mlnation  of  the  dyriamical  matrix  (a  matrix  of  self -consis- 
tently averaged  second  derivatnves  in  t  lie  SCK  theory  as  opposed  to  equili- 
brium vclue.s  in  plain  harmonic  theory)  . 

The  quantum  parameter  A   introduced  in  Appendix  1  is  a  measure  of  the 
quimtum  ir,er},ani  cal  character  (or  zero  point  motion.)  of  a  syste;;;.  Of  the 
tvo  rare  gas  solid.-:  considered  here,  the  value  of  A   for  argon  ±r-    fairly 
small,  vliile  that  for  neon  indicates  somew'iat.  riore  substantial  zero  point 
iiLOti  on .~ "    iherefore  the  uiroerturbed  (i.e.  f:diabatic)  phonon  spectra  of 
argon   and  neon  have  been  ca.1  culatsid  by  the  SC.'i  theory. 

As  is  \'eli  knov-M,  t^iC  p.f'iiicip.'.il  difficulty  i.n  solving  the  adiabatic 
eleci'ropic  ■^truci'vir-^'  p'"obJi='n'.  (7  —  3)  r!''"-'.se.~  froni  t  h"  m:>ny— bod)'  r'P'"';?."e  o^ 
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the  electron-electron  interaction  ccii.bii  e.i  vlth  the  requirement  that  the 
electronic  wave  functions  be  properly  anti-sy-itinietrized.  Almost  all  treat- 
ments have  as  their  aim  the  construction  of  a  ona— electron  model  in  which 
the  electrons  are  treated  as  independent  particles  moving  in  some  well- 
defined  average  field  due  to  all  the  other  electroiia.  Even   at  this  level, 
rigorous  treatment  of  Fermi  statistics  introduces  complicated  exchange 
contributions  to  the  effective  Hamiltonian.  The  inclusion  of  spatial  cor- 
relation introduces  further  complications.  The  net  effect  is  to  make  the 
problem  alniost  intractable  for  a  crystal.  For  example,  correlation  and 
exchange  effects  m.ay  be  treated  through  a  comiputational  realization  of  the 
configuration  interaction  scheme  for  atoms  and  molecules,  but  in  solids 
such  an  approach  is,  to  this  date,  well  beyond  the  limits  of  feasibility 
and  other  schemes  must  be  invoked.  "  A  much  explored  alternative  is  to 
make  soifiC  sacrifice  in  rigor  (small,  one  hopes)  and  thereby  obtain  a  sol- 
uble m.odel  of  electron  behavior.  The  next  section  deals  with  the  form  of 
such  a  m.oG';l  and  the  Hamiltonian  H  (R)  of  (2.-3)  that  results,  with  particul 
emiphasis  or;  the  features  needed  to  achieve  the  eventual  goal  of  tliis  work. 

2-2  .  JI'  he  Xa  I'amilt  o  n  ian 

A  brief  outline  of  the  Xa  method   for  the  non-spin  polarized  case 
will  be  given  here.  Tlie  method  itself   is  v.'cll  documented  and  has  had 
many  useful  applications  in  a  variety  of  circumstances .' ~  For  tb^e  present 
purposes,  'Aie   method   allows  a  coinputationalJ.y  manageable  form  for  the 
Hamil  t:on.iaa  in  calculations  of  the  electionlc  p'ropertief  of  a  solid. 

For  a  static  lattice,  in  tfie  adiebatic  sense,  the  Xu  m.ethod  assumes 
for  tlie  total  electronic  energy  the  expression 


:  >  -  <T,,_>  -1  <:'  .-  (2-c) 

Xa         A-,1         XU: 


where. 


<-^  ^   ^   -   In.  iu^(i)\Cu.(l)Qi.-.  (2-9) 

Xa      .ill    11    —1 


<K   >  =  r'Z."  /R.   -  Vlp(l)(.2Z./r,  Jdr^  +  |  c  (1)  (p(2) /r.  ^)dr^  „ 
Xa     V   J  m  JH'    .■'J        J   Ij   — 1    J  i'^  ~12 

+  1/2  Jp(l)U^^(l)  dr^  (2-10) 

with  Z.  denotir.fc  ths  j '" ''  nuclear  charge  and  R.   the  j    and  m  '  nuclear 
J        -  ]m 

separation.  The  charge  density  is 

p(l)  =  I  n.u*(l)u.(l)  (2-11) 

.11     1 

with  n   the  occupancy  of  the  spin  orbital  u..  The  local  exchange-correla- 
i  1 

tion  operator  in  the  non-spin  poiarl"ed  case  is  defined  as 

U^^/1)  -  -9a(3/8Tr  p(i))^^^  (2-12) 

with  u,   a  paraineter  to  be  cai.ibrated  against. ^  for  example,  atomic  properties. 
Application  of  the  variauloaal  principle  to  (2-8)   leads  to  the  effective 
one-electron   Schr6'd:inger  equation 

[-v:  +  yz./r,.  +  i(2Z_/r,.Jp(2)dr„  -  6a(3/8-  p(l)  )^-^^]u .  (1)  =-  £  u(l)  (2-13) 
1    .'  1   Ij    J    i   i^      -V.  1       11 

The      ouc^-elec  rro'    ei^-envalnes   E.    of    (2-13)    are   related    to    the   occupancies 

1 

n   so  that  Ferni  statistics  are   obeyed. 
1 

The  static  Jattice  cohesive  energy  per  particle  for  a  solid  of  N 
constituent  nuclear  particle,-^  and  cf  a  specified  density  is  defined  as 

(E/N)  =  <£,.  >   ,.,/N  -  <E_>  ^  (2-14) 

Xo:  Eolad       Xa  atom 

This  definitior  is  con£i.-.>tent  with  AppendJ:-.  I  provided  the  two  Xa  total 
encr;;;ier  -ivc    calculated  p.t    the  sa^.>e  '.' . 
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The  selection  of  a  value  for  a  has  been  the  subject  of  much  investi- 
gation and  discussion.  However,  as  a  result  of  previous  calculations 
done  on  argon^  ^  and  neon,^'*  a  value  of  2/3  has  been  used  for  neon  and  a 
value  of  0.72131  for  argon.  The  neon  value  is  that  due  to  Caspar,  Kohn, 
and  Sham,^^  while  the  argon  value  is  such  that  the  virial  theorem  is 
satisfied  for  the   energies  determined  from  a  single  determinant  of  argon 
Xa  orbitals  using  that  v-?lue  cf  a. 

2-3.  The  APW  Method 

The  APW  method  is  a  v;ell-established  scheme  for  solving  a  given  one- 
electron  eigenvalue  problem  for  a  crystal.  An  excellent  review  by  Mattheiss 
et  al.'  discusses  the  self-consistent  field  symmetrized  form  of  the  AF'U' 
method  (SAPW) .  Just  this  form  has  been  used  in  the  present  research. 
Therefore  the  presentation  here  will  follow  the  notation  and  style  of  that 
review.  As  the  calculation  was  done  on  two  solids  vjhich  hax^e  the  same 
sjTTjnorphic  space  group  (see  Appendix  1),  the  explicit  reference   to 
non-primitive  lat.tice  translations  v;ill  be  dropped. 

The  potential,  of  the  functional  form  discussed  in  section  (2-2),  is 
assumed  to  he  spier: cally  symmetric  inside  spheres  centered  at  each 
nuclear  site  in  a    solid  and  constant  outside  thnr.e  spheres.   This  periodic 
muffin-tin  (MI)  potential  is  a  physically  plausible  one  which  approximates 
the  exact  periodic  potential  wliile  being  much  m.orc  convenient.  In'  the 
constant  potential  region  j  the  AI';-.'  basis  functions  are  taken  to  be  plane 
waves;  while  inside  the  sphere  thic  APW  basis  functions  are  expanded  in 
spl.erical  harmonics  and  energy  dependetit  radial  functions. 

The  electronic  Bloch  fui'iction  for  the  j  '   band  and  h  "  wave  vector 
wi''  c-""r<-"  H  ''!'''  m-'v  be  c---'iand'^d  i"  tar-::;  cf  these  AP'.'  basis  functions 
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^.,  (r)    =    y   C(K)$rk+;c,r,E,)  (2-15) 

Jk-         I       -       J 

with   the    sum  K  over   all   recij^rocal   lattxce  vcczors.    The   APW  basis   function 
of    (2-15)    is,    letting  R     be   the  MT   sphere   radius. 


•I)(k+K,r,E.)    -   e:a;[i(k+K)-xi  r>R        (2-iea) 


^j_^(|k+K!R^) 


=   I  1  4Tri  --^xf  IJVk' W^    -£^^'^l)^£j-^)  ^<^s      ^^-^^^^ 

£  m  X-      s      J 

In  (?.-16b),  Y^   is  a  spherica].  harmonic  (and  is  of  the  Condon  and  Sbcrtley 

iCin 

convention"'),  j  „  is  a  spherical  Bessel  function,  P,  ,t^  J^Tid  ^,j_^   a^e,  the 

angular  coordinates  of  the  vector  k+K,  E^  is  the,  energy  eigenvalue  of  the 

i    band  and  k   v/ave  vector  and  the  u,,  (r,E,)  are  solutions  to  the  radial 

1^  J 

SchrJidinger   equation    (in  R^^dberg   units) 

-1/r^   d/dr    (r^   du^^/dr)    +    [  (£(£+1) /r'')    -r  V^,^J    u_^    ---   E.Uj,    (2-17) 

where  V   is  the  muffin-tin  potential. 

The  Bloch  function  of  the  j    band  and  k   wave  vector  can  be  syFivoe- 

th  th 

trized  to  transform  irreducibly  as  the  n   basis  partner  of  the  a   repre- 
sentation of  the  group  of  k.  First,  projection  operators  are  applied  to 
the  Al'W  basis  functions  of  (2-16)  to  obtain  the  the  SAPW  basis  functions 


n,  s 


(k4-K,r,E")  -  I      r*"^(R)Ri(k4K,r,E.)  (2-18) 


ReG, 


.    -     r  •  ^        th         ,      th      ,  ^      ,         th        ^     .  ,       .  .   _, 

v;itn   1^'         the   u        ana    s        exfc.meiit    of    the   a        matrix   representatJ  on   or    R, 
n ,  K 

and   R   belonging    to    the    operations   of    the    group    of   k    (G    )  .    Then,    with    t?ief 
SA}"/,    the    syim^.ecrized    Blocli    functions   are   obtained 

^^^(O    =   11^;''      (R)v''      (h-K,r.E'b  (2-19) 

nk  —  f;    ^_      n ,  s   -  -     n ,  s   — 

n:  th 

with    E'    "in    hotVi    ("^-18")    anfi    (')-iQ')    *.he    eneipv    eigenvalue   of    the   O. 
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represent&rion .  Strictly  speaking,  thi-S  eigenvalue  should  have  a  basis 
partner  index,  but  for  tr,ost  cases  (as  in  the  present  calculation)  the 
eigenvalue  is  degenerate  in  the  basis  partner  index. 

The  effect  of  a  rotation  operation  R  on  the  APW  basis  function  (2-16) 
is  defined  by 

R  *  (k ,  r ,  E  .  )  =  0  ( k ,  R ""'■r ,  E .  ) 

=  $(Rk,r,E^)  (2-20) 

with  R        the   inverse   operation   to  ?~_.    The   final   form  of   the   SAPW  becomes 

*'^      (k-:-K.r.E^)    =      y      r*^    (R)    exp[lRk-r]  r>R      (2-21) 

n,s REG,    ^'"'"■"  ~^~ 

oumj  -  „ 

=      y      r*^     (R)    y   y    Airi"-- —     Y^:    (ei   ,  ,o  ,  )    u„(E",r)Y,     (e,ct)    r<R      (2-22) 

Rr.G,       '  i  m  u,  (E    ,R  )  —     — 

—     k  is 

where  on  the  R£S  of  both  (2-21)  and  (2-22) 

Rk  =  R(k+K)  and  k  = [Rk |  =  iR(k+K) |  (2-23) 


In  most  applications  using  the  APW  basis,  the  addition  theorem  of  the 
spherical  liarnionics  (A3-12)  may  be  invoked  to  further  simplify  (2-22) 

^"^      (k4-K,r,E")  =  A7r  V  r*°  (R) ^  i^  -^^-^'-^  P.  (=3^)  u„  (E",r)  (2-24) 

^'^"-'-  Ren  ^-^~£    u„(E^.E  )  ^^  ^^     ^- 

-~  k  £     s 

The  purpose  of  this  syTiur.etrication  is  two-fold:  the  reduction  of  the 
size  of  the  secular  matrix  to  be  solved,  and  clear  identification  of  the 
coni:ribut:ir;:\'.  SyTrrr.etry  to  the  electronic  bands.  As  a  consequence  of  using 
SAPV.'s  ,  the  entire  problem  m.ust  be  foi-mulated  in  symmetrised  terms.  This 
means  that  the  phonons  of  the  nuclear  m.otion  for  the  clectron-phonon 
interaction  considered  here  must  be  expressed  in  terms  of  the  appropriate 
irreducible  representations  and  basis  partner  (which  is  equivalent  to  the 
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phonon  wave  vector  and  polarization)  .  The  results  cf  the  Sj-'mT!ietri?.ation 
of  the  phonons  (in  the  form  of  their  dispersion  ciarves)  are  given  in 
Figure  2-1  for  the  typical  fee  solid  with  one  atom  per  unit  cell.  In 
Table  2-1  are  given  the  eigenvectors  of  the  phonons  at  each  of  the  sirnnne- 
try  points  and  at  each  polarization  (or  basis  partner) . 

One  final  point  of  the  A?W  method  concerns  normalization.  Pioper 
normalization  is  required  if  the  elt'ctronic  Bj  och  iunt.tions  are  to  be 
used  outside  of  the  actual  SAPW  program.  The  normalization  integral  of 
(2-19)  is,  in  general,  over  the  entire  periodic  macro-ceil;  however,  using 
Bloch's  theorem  the  integral  may  be  reduced  t5  N  times  the  integral  ever 
a  single  unit  cell  (where  N  is  the  number  of  unit  cells  in  the  periodic 
volume).  Reference  3  shov7s  that  the  application  of  standard  group  tlieore- 
tical  methods  gives  for  the  normalization 

'   —      —  K   K  K£(j 

j  (!pi_k'-(k+K)lR  )         i;(k'R^)  Jo(kR  )  ('^'s  p 

4^p/  ~^-l^- ^:^ -^-  +  47Ty(2a+l)-^:: ^     u*u,r  dr]  (2-25) 

I —   i  g        I  usCK  ,K  ;n,,  (1  ,R  )    o 

k  s   X      E 

where  g  is  the  order  G,  ,  d   is  the  dimensionality  of  the  a.        representa- 
°  k   a 

tion,  and  fi  is  the  primitive  unit  cell  volumt;.  The  notation  of  (2-23)  has 
been  used  in  (2-25)  and  in  Appendix  III  where  the  formulae  tor  the  EPME 
are  derived  based  on  the  SAPW  in  the  form.  (2-22)  . 

2-A.  The  Fr^hlich  Kr-mi]  tonian 

Though  Frtihlich"*' '^  and  many  others'"^  using  his  approach  had  the 
principal  objective  of  explaining  m.etallic  superconductiv'ity ,  the  Hamil- 
tonian  and  approach  to  the  electron-phonon  interaction  which  he  form.ulated 
is  of  a  Reneral  nature  (within.  oF  course^  certain  approxirnations)  and 
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Table  2-1 


Phonon  Eigenvectors  of  the  Six  Non-equivalent  BZ  Points 


SycTiaetry  Matrix  Form  of 

Point  Eigenvector 

L       Tl     T2 

1 

r  0 

0 
0 

A  i 

0 

0 
X  1 

0 

li 
I  li 

0 

1      1  -  li   1  +  li 
/3     'i     A2      2      v^2 

/3      2   /r2   2'  /i2 
1        li     -li 

T       [^  R 

v3       v^i      v-i 

0       0       i 

V:  2    2    2    2    0 

1  +  li  -1  -;-  li 

2  1        2        2  0 

^'   Row  dtinotep  x,  y,  z  conponeats:  column  deuoLes  polarization 
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The  first  two  eigenvectors  of  W  are  actually  pseudo-longitudinal 
and  the  ].ast  is  actually  pseudo-transverse. 
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will  be  utilized  in  this  work.  Ari  excellent  paper  by  Hedin  and  Lundqvist 
(the  author  is  grateful  to  Dr.  N,Y.  Ohm  for  bringing  this  reference  to  his 
attention)  gives  a  strictly  formal  justification  for  the  assumed  fonn  of 
the  FrHhlich  Haniiltonian  that  has  been  used  here.  The  purpose  here  is  to 
present  that  assumed  form  within  certain  approximations  and  shov;  the 
natural  connection  to  the  Born-Oppenht;inier  Hamiltonian. 

As  a  convenient  starting  point,  the  Hamiltoriian  H  of  (2-2),  for  N 
nuclei  of  mass  M,  is  expanded  about  the  equi.libriuiri  lattice  sites 
(R   )  in  powers  of  thfr.  displacements  (5R^)  .  In  keeping  with  the  hamionic 
model,  the;  gradient  of  the  total  energies  of  the  electrons  and  of  the 
nuclei  taken  togeclier  vanishes.  The  only  tertas  left  to  first  order  in  the 
displacem.ents  are  the  electron-electron  ,?nd  nuclear-electron  interactions 

H(R°)  =   H  (P°)  +  H  (R°)  +        H^^„ 

=  [T  (R°)+V  (P°)]  i-  [T  (R^')+V  (r")+V   (R°)1  +  5R-V/  (V  +-V  )  I  o    (2-26) 
n  —    n  —        e  --    e  —    ne  —  '      —  ix  ne   e  _R 

H   is  the  electron-p'nonon  interaction  to  first  order  in  the  displacecient 
and  by  an  obvious  simplification  of  (2-2),  T  has  replaced  the  kinetic 
energy  suras.  The  individual  operators  of  (2-26)  ir.ay  be  expressed  (see, 
for  instance,  references  22  and  31)  in  second  quantized  field  operator 
form  as 


H  (K°)  -    y   E^a'\a  . 

-   I      (-,.  +  1/^)  <^  (2-27) 

51?'   " ' 


with  the   a's   denoting  boson  quasi-particle  operators 
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--=  J  u,  E°  (2-28) 


with  the   c's   denoting  fermi  quasi-particle  operators,  and 

H^^=  J^^M(k';k)c^,c^(a:^+a^)  (2-25) 

In  (2-29),  the  sum  is  restricted  to  ^  =■•  k'  -  k  +  K  (K,  a  reciprocal 
lattice  vector),  the  displacements  have  been  expressed  in  normal  mode 
coordinatsc 
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6R.  -  (NM)  ""-^^  I      Q  ,exp[i£-R"]c  , 

=  (NM)-^/^  I^  (20^^)-^/^  (a:^.-a^)        (2-30) 

and  the  electron-phonon  ratrix  element  (EPME)  [between  electronic  B.loch 
states  (2-3)]  is  defined  as 


M(k';k)  =       <V'^-EP'V 

=    (N>;2to'',)~^''^<^'ME   a'^p^V     +V)j,oi>Ff>         (2-31) 
c^k  k       <^'.      R     ne     e     ii       iL 

The   energies   of   equations    (2-27)    and    (2-28)    have   the   same  meaning  as 
those   of    (2-3)    and    (2-7)    and    the   superscript;      o     means   the   static   equil.i- 
briuTii   lattice    (unperturbed). 

The    totality  of    (2-27),     (2-28),    aaa    (2-29)    comprise   the   assumed    form 
for    the   Fr'Jhlicii   Harailtonian  K     used   in   this  vrork 

H      =   ;    E%;'c.    +  J      E^a\a    .    +     ^     M(k' ;k)c;^ ,  c.  (a"^     +a      )  (2-32) 

F        i-      k   k    K        '',       c,l   qA    cja        /^    ,      --     -     k      k      -3A      ^A 
k      — qA     -      ^i      -L  k,^  " 

Some  explanation  cf  these  equations  is  no\j   in  order.  In  (2-28)  and 
(2-29),  the  k  and  k'  also  ir.ply  baud  indices.  Bra  and  ket  notation  simpli- 
fy the  electronic  coordinate  intej^rration  in  (2-31).  The  Bloch  condition 
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of  (2-1)  has  been  used  in  (2-31)  to  remove  trie  exponential  of  (2-30)  and 
to  give  the  selection  rule 

1  =  k'  -  k  4-  K  (2-33) 

which  governs  the  sum  restrictions  in  (2-29)  and  (2-32).  The  presence  of 
the  reciprocal  vector  K.  allows  for  uraklapp  processes.  ^'^^  The  convention 
of  the  wave  vectors  throughout  this  work  is  based  on  (2-33)  with  k_  denoting 
the  initial  state  (valence  or  "p-like"  band),  k'  the  final  state  (conduc- 
tion or  "s-like"  band),  and  q   the  phonon  wave  vector. 

The  eigenvectors  and  normal  mode  am.plitudes  cf  (2-30)  are  chosen  to 
have  the  foJ.lowing  phase  convention 

Q^A  --   '^-s^X  (2-34a) 

r:  ,  =  c   -  (2-34b) 

The  eigenvectors  for  the  syirjuetry  points  used  in  this  work  (see  Appendix  I) 
were  given  in  the  previous  section  (Table  2-1)  and  have  been  a.dapted  in 
symmetrized  form  to  the  present  calculation  from  those  vectors  of  refer- 
ence 32.  These  eigenvectors  ate  seen  to  satisfy  the  following  orthogo- 
nality and  closure  relations 

e   -E   .-  0  (2-35a) 

— (.^A  "~C[A       A  A 

I  t:'\  c^,  =  6.  .  (2-35b) 

with.    i,j    the    X,    y,    or    z    conipcaent    of    the    eigenvector. 

In  writing   down   11,   of    (A'-?>2)  ,    the   explicit-    electron-electron   and 
e]  ectrc)n-nuc].ear    interactions   are    ai.:counted    for   by   using   an   appropriate 
screer.ed    forr:   of    the   poteni.i,'ls    in    tlie   EPKE   expiession    (2-31)    and   by    cal- 
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dilating  the  "bare"  or  unperturbed  phonons  cf  the  Harailtonian  R^  using 
some  method  that  gi\'es  the  correct  small  wavenumber  behavior  of  the  phonon 
energies  (ji.e_.  proportional  to  the  wavenuinber)  .  The  screening  is  inherent 
in  the  Xa  potential  approximation  for  exchange  and  correlation  effects 
and  the  correct  behavior  of  the  phonon  energies  is  contained  in  the  SCH 
theory  if  any  reasonable  (though  not  necessarily  accurate)  total  lattice 
potential  (,in  this  work,  the  Lennard-Jones  superposition)  is  used  to 
calculate  the  unperturbed  phonon  spectra.  Therefore,  as  used  here,  the 
bose  and  ferir.i  creation  and  annihilation  operators  of  (2-27)  and  (2-28) 
are  formally  "quasi-particle"  operators  as  opposed  to  the  usual  single- 
particle  electron  or  phonon  operators.  This  peine  is  emphasized  and  given 
a  strictly  formal  justification  by  Fedin  and  Lundqvist.^°  H,,,  in  the  fore 
of  (2-32),  therefore  will  be  assumed  to  be  applicable  to  the  present  cal- 
culation. It  should  be  mentioned,  however ^  that  energy  corrections  cal- 
culated using  H^,  have  been  found  to  affect  only  electrons  cJose  to  the 
Fermi  surface,  v.'hlle  the  phonon  dispersion  v/as  altered  substantially .  ^" '  ^^ 
These  facts  vjill  be  used  both  as  an  aid  and  as  a  check  for  this  calculation 

By  Eiea;:is  of  H.„,,  the  energy  corrections  may  be  calculated  using  second 
order  Raylcigh-Sciircdluger  perturbation  theory.  First,  we  define  some 
quanticles  for  purposes  of  simplification 

H^  =    R  (R°)   +    H  (R^) 
u  -  e  — 

=  J  E^-'a'a    -i-  y  E^'c'c,  (2-36) 

5.  h.     

o    r   o         V      n  '  . 

E  =  1    \   n    +  J  K^,  n^^  (2-37) 

wherc'  the  polari/.atior.  labf-1  is  Implicit  in  ^  as  is  the  band  index  in  ]:_, 
aiid  u,  and  v.      d^noze    the  occuBation  nu;:.bfcr:_s  of  the  electronic  state  4c  and 
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phonon  state  c^_   respectively.  It  should  be  noted  that  the  electronic  occu- 
pation numbers  have  an  implicit  spin  suitjination  factor  of  twc  contained  in 
thein  and  that  both  the  phonon  and  electronic  occupation  numbers  are  not 
necessarily  integers  because  of  the  "quasi-particle"  nature  of  the  field 
operators.  The  total  energy  to  second  order  in  the  perturbation  H_  Tr^y 
now  be  ^^^ritt^n  as 

E  =  E°  +  <H|Hgp|=>  +  <H|H^p(E°  -  H°)-\p!H>  (2-38) 

where  bra  arid  ket  notation  indicate  nuclear  and  electronic  integration 
and  E   is  the  ground  state  total  systen;  wave  function  [in  the  sense  of 
section  2-1  and  (2--4)]  with  n^  electronr^  in  state  k  and  n  phonons  in 
state  q^.  The  first  order  energy  terra  in  (2-38)  vanishes  due  to  the  haraonic 
nature  of  the  phonons  (foriBally  this  is  equivalent  to  the  creation  or 
destruction  of  a  plionon  in  the  tota].  syste~  and  the  resulting  ortho:?o- 
nalitj'  of  that  system  state  to  the  systen  ground  state) .  The  total  second 
order  energy  correction  of  (2-3S)  may  be  shown  to  be" 

E*^^^  =  <H;|  y   |M(k';k)rra''  c;'c.  ,(E°-.H")"^a   c.^.c, 
y"^],i   '   -  -  '    -q  k  k  -a.  £  ll 

+  aVc.,(K°-K^)-^a  cj  ,c,  ]\->  (2-3S) 

q^  Ic  k  ^  ^     ^ 

where  the  sum  is  restricted  by  (2-33).  For  clarity  and  illustration,  Figure 

2-2  shows  the  intieractions  inherent  in  the  assumed  fonri  of  F.  ^  and  the 

tr 

processes  that  contribute  to  tne  total  second  order  energy  of  (2-39).  The 

use  of  (2-36)  and  (2-37)  in  ter!:is  of  occupation  nu":b:'rs  reduces  (2-39)  to 

n_^  (n  +1) 

F^2)  .  Z   |K(}^';il)IX'a-.,)[^-~-;;-H-  -:~^-~--l  (2-40) 

k,k'  -     -  E,  -i'.,  ,-H,:       E,  -E  ,-u 

where,    again,    the   suju   is    restricted   b}     (2-33).    liie    synunetry   of    the    ei.gen- 

veciors    ("L^-HAbl    <"Tves  (■■}     =   ''      :    t?",;^    in    tl'^r-.'J    P3ull''hr:  l!"i,    it    is   cxpec- 
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Figure  2-2.  II]  ut.trafioa  of  liie  interactions  involved  in  th.j  perturbation 
Harciltonian  i!  _  r.nd  of  the  resultant  second  oi'der  energy  pro- 
cesses.   The  various   pictures   and    their   correspondin.g    teinns    are 
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ted  that  n  =  n   .  Therefore,  (2-40)  may  be  simplified  to 
'1-3.  '.        ^ 


(2-41)  is  the  total  second  order  energy  correction  to  the  system  and  the 
electron-phonon  effect  on  the  phonons  is  in  the  tern;  proportional  to  the 
phonon  occupation  n  . 

Identification  of  the  phonon  perturbation  with  the  total  energy  E  of 
(2-38)  necessary  to  increase  the  q^   '  phonon  occupation  by  n   gives^^ 

=  3E 
qX     9n  , 

with  the  sutti  over  Ic'  still  restricted  by  (2-33)  and  the  polarization 
explicitly  included  as  a  reminder  that  the  transitions  are  between 
representations  of  the  wave  vectors.  For  this  calculation,  the  sum  over  t' 
is  only  over  available  scattered  states  in  the  conduction  or  excited  bands, 
Elxpressiou  (2-42)  has  been  programmed  (see  Appendix  IV)  to  calculate  the 
total  sccoiu;  order  energy  corrections  only  for  transitions  involving  all 
"p-like"  filled  valence  b;md  states  to  all  excited  empty  "s-lJ.ke''  conduc- 
tion band  states.  These  ti'ansitions  are  based  on  the  selection  rules 
calcu].ated  in  Appendix  IT  between  the  electron  wave  vector  representatioviS 
and  the  phonon  v?ave  vector  representations. 
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2-5.  Tne  Derivation  of  the  'Forj.  for  the  EPME 

In  section  (2-4),  the  existence  of  an  EFME  was  established  and  shoxvn 
to  follow  as  a  natural  consequence  cf  the  l^om-Oppenheiiner  theory  and  was 
used  in  the  Fr8nlich  Hamiltonian  to  calculate  second  order  energy  correc- 
tions to  the  phonons.  The  purpose  of  this  section  is  to  derive  and  justify 
an  e:Mpliclt  form  of  the  EPME  specific  to  the  APt\'  formalism. 

The  muffin-tin  (MT)  potentials  are  assumed  to  move  rigidly  vrith  the 
nuclei.  This  assumption  is  similar  to  the  "rigid  ion"  approximation, 
the  fundamental,  and  physically  more  reasonable,  distinction  being  that 
the  MT  potentials  are  a  self-consistent  description  of  the  Rorn-Oppenheimer 
solid  rather  than  a  simple  superposition  of  atomic  potentials.  Because  the 
MT  spheres  are  traditionally  chosen  to  touch  so  as  to  maximize  the  frac- 
tion of  the  cell  volume  they  enclose,  there  appears  to  be  a  problem  of 
overlapping  and  mixing  of  adjacent  potentials  when  the  Mils  are  displaced 
(a  point  originally  em.phasized  to  the  author^"'').  Kenney^''  has  stated  that 
the  validity  and  applicability  of  the  APW  method  are  rather  independent 
of  the  sphere  radius  provided  the.  sphere  rem.ains  v/ithin  che  Wigner-Seitz 
cell  (primitive  unit  cell).  Fe  also  gives  a  plausible  argument,  in  the 
case  of  mietals,  that  the  sphere  radius  be  reduced  five  to  ten  percent  from 
half  the  nearest  neighbor  distance.  Ziman^^  devotes  a  considerable  dis- 
cussion to  the  effects  of  large  displaceinerits  (sm.all  wavenumber)  which 
is  really  only  pertinent  to  metals  or  ser.i-conductors.  Estimates,  based 
on  Keuney's  suggestiou,  of  the  mean  square  displacement  (in  Rydberg  atom.ic 
un  its) 

(•■.Y   -  1  /(M:^^-^  (2-^3) 
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for  neon  and  argon  show  thsr  the  average  displscements  are  about  fifteen 

percent  for  neon  and  about  ten  percent  for  argon.  In  spite  of  this,  and 

due  to  the  fact  that  Z±m£:n   and  Keaney  consider  large  displacements  only 

for  TTiGtals,  the  overlap  effect  has  been  neglected.  Both  the  smallness  and 

the  flatness  (small  derivative)  of  the  potential  in  the  overlap  region 

(as  contrasted  to  metals)  tend  to  make  this  approxicuition  plausible. 

Following  Golibersuch, ^  the  Bloch  function  for  the  j    band  at  the 

k   wave  vector  in  the  perturbed  lattice  (i_.e_-  perturbed  by  the  electron- 

phonon  interaction)  is  given  approximately  by  first  order  perturbation 

theory  as 

M(k';k) 

m^    ^  T°  +  y t"  (2-44) 

-    --   k'  E°-F°   -^ 
-   k  k' 

with  the  band  indices  implicit  (as  in  section  2-4)  in  tb^e  wave  vector  and 
wi.th  th.e.  suR.v.iation  in  general  over  all  states  including  core  and  valence 
(occupied)  and  conduction  and  excited  (unoccupied) .  The  superscript   o 
vrili  denote  throughout  this  section  the  unperturbed  static  equilibrium 
lattice,  v/hlle  its  absence  v/ill  denote  the  perturbed  lattice.  Th^e  explicit 
_r  dependence  and  the  band  index  j   of  section  2-3  will  be  dropped 
henceforth  (except  where  essential  for  cjar::ty)  with  the  convention  that 
the  wave,   vector  Includes  the  band  labeling.  The  square  of  M(k';k)  is  for- 
mally the  first  order  (in  the  wave  function)  transition  probability  be- 
tween the  eqi'.illbrium  states  k  and  1^'  v:here,  as  in  section  2-4,  Ic  is  the 
initial  state  and  k_'  the  final  state.  It  is  therefore  appropriate  to  take 
this  quantity  as  the  ET'ME.  Froc;  the  viewpoint  of  scattering,  tVie  EPME  can 
be  related  to  that  part  of  the  transition  matrix  (  T   matrix)  linear  in 
the  nuclear  displacement. 

The  static  nTippr'"\irbed  RToch  * ''ncti  on'^' ,  simplified  P""d  superscri  veQ 
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from  (2-15),  become 

"^1  =   1   C°(k,k^K)-:^"(krK)  (2-45) 

k   ^ 

where  the  coefficients  now  contair?  aa  explicit  reference  to  the  initial 
state  wave  vector  k.  This  exjxlicit  reference  is  to  remind  us  that  the  per- 
turbed wave  function,  constructed  frnr.'.  "-.he  ]-ertuited  APW  basis  functions, 
is  not  necessarily  a  periodic  function 

\  "  I   C(k,k-!-K)0(k+K)  ^    y    C(k,k')$(k')         (2-46) 
-   K    '        '  '  k'?ll+K 

where  the  k  label  on  the  LHS  is  only  to  indicate  the  initial  eigenstate 
of  the  systeiTi. 

In  (2-46),  all  possible  wave  vectors  are  needed  due  to  the  non-peri- 
odic nature  of  the  perturbed  wave  function;  that  is,  the  label  k  no  longer 
denotes  eigenstates  of  the  now  perturbed  system.  Expanding  the  perturbed 
APV7  basis  functions  and  coefficients  of  (2-46)  to  first  order  in  the 
nuclear  displaceinents  gives 

C;(k,k+K)  =  c''(k,k-:-K)  +  CR-V'  c!  (2-47) 

— ' — ' —   R   -r.O 

—   K 

^(k+K)  =  ^"(k+K)  +  £R-V^(i>|  ^  (2-48) 

where  6R  are  the  displscexeuts,  V   are  the  gradients  v.'itb.  respect  to  the 

nuclear  coordinates,  and  F^  are  the  equilibrlurt  configuration. 

Using  (2-47)  and  (2-48)  in  (2-46),  and  keeping  only  constant  and 
linear  terms 

\  -  ):c°(k,k^)0^h+K)  +  yc"(k,k«;)V^c!  ^  +  [V^^C!  ^""(k+T)   +   VVj^d  ^<°(L') 
-   K  K  "-  R    K-'  Ji  kVk+K" 

-  •<  +    J        V  Cl  .0"(k')  +  y  C°(k,k-IK)\\,d!  ^  (2-49) 

-   alT  k'   -  R*"         K  "  R^ 


The  second  terci  of  (2-49)  represents  scattering  into  an  available 
equilibrium  eigenstate,  since  the  sur.i  i?  ever  all  possible  k'  states  and 
reciprocal  lattice  vectors  K.  "Availabl?"  is  emphasized  because  all  scat- 
tering is  considered  from  the  ground  state  configuration,  and  the  lowest 
energy  eigenstates  are  filled  to  the  Fermri  level  and  are  NOT  available  as 
final  states  into  which  scaLterin:^  occurs.  Rewriting  the  second  term  to 
emphasize  this  fact 

--   all  k'  -  R^ 

=   y   A(k';k;+°  (2-30) 

;  ,    —  '  -  k 
K_         — 

where  the  sum  k'  over  the  scattered  states  is  now  restricted  to  the 

excited  or  conduction  band  eigenstates. 

The  third  tern;  of  (2-49)  is  seen  [through  the  unperturbed  coefficients 

C°(k,k^K)]  to  allow  the  unperturbed  APW  basis  function  to  follow  the  motion 

th  th 

of  the  nuclei  and  the  rest  (those  not  in  the  j    band  and  k   vector  eigen- 

state)  of  the  electrons.  This  is  also  pointed  out  by  Kenney^   in  a  dif- 
ferent context.  For  convenience,  this  term  will  be  rewritten  as 

B  =  I  C°(k,k+K)V  e]  (2-51) 

--   K  --  R 

Comparing  (2-44)  and  (2-49)  shows  that 

M(V'  •  !,-> 

——"    -   <^F°  lE  >  +  <-°  Is  >  (2-52) 

^o  _o      <*-   ii      ii   i:- 

The  electronic  Sc!i:r?3 dinger  equaticn  for  the  perturbed  situation  may 
be  taken   fro-;  (2-3)  as 

H  T,  =■■  F^  V,  (2-53) 

e  k    k  k 
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with  V    (r)    and  V      (r,R)    the.   sum  of   displaced   rnuffin-tin  potentials 
e  —  ne 

V      (r,R)   +  V    (r)    =    I   V_(r-R^)  (2-54: 

3 

ar!.d  j  labeling  the  nuclei  and  their  respective  muffin-tins. 

Expanding  (2-53)  and  (2-54)  about  the  R,  equilibrium  positions  to 
terms  linear  in  the  lattice  diaplacerr.ents  (5R.),  one  finds 


H°(r'')[S,  +  B,  ]  +  6V'f°  -  £,"(3,  +  B,  )  (2-55) 

e  —    k    k       k    k  k    k 


9V 

vnth  H  (R  )i'  =  E  w,   and   5V  =  >  oR.'-^- —  ^. 
e  —   k    k  k  ^  — J  dR.  '  ,o 

(2-50)  and  (2-52)  may  ba  us^d  tc  obtain 

M(k';k) 

=  A(tc';k)  +  <H'°  |b  >  (2-56) 

o   o      -  '-      k"  k 

E  -E  ,  —   — 

£  ii 

Solving  for  A(k' ik)  by  vsing  (2--55)  gives 

(e"  -E°)A(k';k)  =  (E°-E°,;<,!b  >  -  <?°  i6viY°>       (2-57) 


With  the  use  of  (2-57)  and  (2-5G),  the  final  expression  for  the 
EPrlE  becoraes 

M(k-k)  .<7;;JI6R^--^M!  J,°>         (2-58) 

In  Figure  2-3,  the  "rigid  ion"  assumption  is  illustrated,  shov;ing 

that  the  potential  at  any  point  within  the  j    muffin-tin  is  the  same, 

v7ith  lespect  to  the  center  of  that  muffin-tin,  regardless  of  the  displace- 
ment 


r.  =  r^  -  r"  =  r  -  R.  (2-59) 

-^        -^3    "3    -   -^ 

with  5R  =  p   -  R"   and  that  as  a  consequence  the  potential  gradient 

~3   -j   -^.i 
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Equi]  ibriuru 

.th  „  ^^. 

J    Muffin- tin 


Displaced 

.  th 

J   Murran-tm 


Origin 


Figure  2-3.  Denon.sLratioa  of  ri^^ici  ion  approximation  for  nuf fin-tin 
diaplaceiiient.  Tliis  figure  shows  that  the  potential  v.'jth 
respect  to  the  muffin-tin  center  is  the  same  regardless 
of  di.snlacement . 


33 


of  (2-58)  may  be  written  as 

_-iiL-J-     =  -  ^^_J_ ^-±L-ZJ.     (2-60) 

3^j       R?     '^  R°  or"-' 

■^       -J  -^      - 

Since  the  MT  potential  is  spherically  synuuetric  within  the  sphere  radius 
and  constant  outside  the  sphere  radius,  (2-60)  becoir.es 

W^(r)  =  e  4^  for  r<R 

MT  —     r  dr  s 

=  0  for  r>R      (2-61) 

s 

where  e   is  the  unit  radial  vector, 
r 

In  section  2-4,  the  nuclear  displacements  (cR.)  were  expanded  ±n 
normal  mode  coordinates  for  use  in  tha  assumiid  EPME  expression  (2-31). 
Comparison  between  (2-31)  and  (2-5S)  shows  thcJt  the  ^iJ'VJ  fo:cm  of  the  EPME 
may  be  written 

M(k';k)  =  (MN2aj'',)"^^^  |   <,(r)<:   -V^^  V_^°(r)  dr    (2-62) 

j^^  unit  cell  -^  - 

where  G  ,  is  the  eigenvector  and  lo  ,  is  the  energy  or  frequency  of  the 
-3_A  ^A 

q  phonon  of  A  polarization,  H  on  th.^  RHS  is  th^-  nuclear  mass,  and  N  is 
the  number  of  nuclei  cr  unit  cells  (for  a  Bravais  lattice).  Using  (2-61) 
and  (2-60),  the  final  expres:;ion  for  the  E?ME  btccries 

M(k';k)  ^  -(M>;2o^)-^/^  f  ^lAD^x'^^l^^^L        (2-63) 

-^       { tti  --    .^ 
J    sphere 

where  the  selection  rule  of  section  2-4  (2-33)  or  Appendix  II  (A2-3)  raust 
be  satisfied.  Appendix  II  shows  that  this  selection  rule  allows  only  cer- 
taiti  transitions  between  the  representations  of  the  electron  wave  vectors 
and  the  polarizations  of  the  phonon  wave  vector.  In  Appendix  III,  explicit 
formulae  are  derived  for  the  SAH.'  basis  function  (2-22)  form  of  (2-63). 


CHAPTER  III 
RESULTS  OF  CAI.CITLATION 
3-1.  APW  Energy  Band  Calculation  of  Argon  and  Neon 

The  Xa  energy  bands  and  wave  functions  in  argon  and  neon  have  been 
calculated  for  the  two  highest  occupied  valence  bands  and  the  lov^ast  mi- 
occupied  conduction  band.  For  neon,  this  would  correspond  to  tVie  "2s-like'' 
"2p-like",  and  "3s-like"  bands;  and  for  argon,  the  "3s-like",  "3p-like", 
and  "4s-like"  bands.  The  other  lower  lying  core  bands  were  treated  in  the 
usual  fashion  as  states  of  the  crystalline  potential  having  atomic-like 
character;  thus  they  correspond  to  flat  (i.£-  k-independent)  bands. 

As  is  expected  for  an  insulator,  the  valence  bands  of  argon  and  neon 
are  conipletely  filled  with  a  definite  energy  gsn  between  the  highest 
valence  band  and  the  bottom  of  the  conduction  band.  By  construction,  the 
bands  also  exhibit  the  features  characteristic  v)f  fee  structures."" 

All  band  calculations  were  done  at  nineteen  s>-r_raetry  inequivalent  k 
points  in  the  irreducible  wedge  (see  Figure  A2-1  and  Appendix  H) .  This 
corresponds  to  256  possible  wave  vectors  or  a  4x^x4  periodic  macro-cell 
(64  unit  cells  v/ith  4  atoms  par  unit  cell  fc  the  fee  structure).  The 

angular  inomentum  sums  v;ere    evaluated  through  x,=10  and  the  riaximun  square 

'>      2 
magnitude  used  for  a  reciprocal  vector  was  SOr  /a  .  Iteration-to-d.terat ion 

convergence  on  the  total  eneigy  was  typically  at  least  10    F.yd .  and 

usually  better.  As  a  result  of  this  energy  convergence,  the  pressure, 

virial  ratio,  and  •^oten^'ir']  r-lso  v/ere  well  converged  wliicli  ensured 

numerically  the  best  possible  eigenvalues  and  eigenf unctions . 
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Relevant  results  of  the  neon  calculations  are  given  in  Table  3-1. 
In  that  Table,  "X- Alpha  Equilibrium"  refers  to  the  static  lattice  constant 
detemiined  to  correspond  to  the  lowest  crystalline  total  energy.  "Experi- 
mental Equilibrium"  means  the  calculation  done  at  the  experimental 
constant. 

Neon  at  the  calculated  equilibrium  lattice  constant  (7.7  a.u.)  is 
seen  to  be  overbound  (in  comparison  with  experiment)  b}''  a  factor  of  two. 
This  overbinding  is  also  evident  in  the  valence  and  conduction  bands.  The 
valence  band  ^-ridtlts  are  expanded  by  nearly  a  factor  of  two  as  compared  with 
the  other  neon  (8.4348  a.u.)  calculation.  The  bottom  of  the  excited  "3p-like' 
ccvnduction  band  extends  into  the  excited  "3s-like"  conduction  band;  in  all 
the  other  calculations  (including  argon) ,  the  lowest  excited  conduction 
bands  are  well  separated.  This  extension  would  add  in  principle,  another 
representation  to  the  "s-like"  representations  of  Table  A2-1 ,   v^hich  were 
used  in  the  actufil  selection  rules  between  the  "p-like"  valence 
(occupied)  and  the  "s-likc"  conduction  (unoccupied)  bands.  However,  the 
additional  contribution  was  not  considered  in  the  calculation  of  the 
phoncn  energy  corrections  and  only  the  "s-like"  representations  of  Table 
A2-1  were  used  in  the  calctilat-ions. 

The  argon  results  of  Table  3-2  appear  to  be  very  good,  particularly 
the  calculated  cohesive  energy.  However,  the  caJ.culated  band  gap  is  only 
about  slxt^r'  percent  of  the  experinenLal  va.lue. 

Since  tills  vrork  is  not  a  coinparativc  study  of  one-electron  models, 
but  an  attempt  to  extend  a  particularly  familiar  example  of  such  a  model 
into  a  previously  unexplored  area,  a  discussion  of  the  relative  merits 
of  the  Xa  calculations  is  not  necessary.  Rather,  the  best  possible  solu- 
tions of  the  X'i-Liodel  that  can  be  obtained  vrill  be  utilized,  regardless 
of  experir.iei'tal  cor^.parison  . 
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Table  3-1 


Calculated  and  Observed  Properties  of  Neon 

Calculated  Observed 

X-Alpha       Experimental  Experimental 

Equilibrium      Equilibrium  Value 

Lattice  Constant              7.7           8.4348  8.4323 

Atonic  Mass   -,oNe^'^       (3.644308xlo'*)   (3,644308xlo'^)  3.644308x10  ' 

Exchange  Parameter-Alpha       2/3             2/3  - 

Pressure(in  kilobars)          -0.24          -4.59  0.0 

Total  Energy  -254.9852860    -254.9844249 

Total  Atomic  Energy  ^     -254.9815919    -254.9315919 

Cohesive  Energy            -0.0036941      -0.002833  -0.0019186 

Zero  Point  Energy            9.3644>:10~^    7.5127xlO~^  4.90S2xlO~ 


Lennard- Jones  Parameters 


-4   .  „__  _-4     .  .., -4 


e  4.29037x10     3.29028x10       2.2990x10 

a  4.9S436556     5.47097073       5.2695 

Quantum  Parameter  '^  5.128128x10"^   5.572535x10"-^    8.3457xl0' 

^  All  quantities  are  in  Rydberg  atomic  units  except  where  no'.ed 

Gee  references  21  and  36 

See  reference  24 

See  Appendix  1  for  definition 
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Table  3- 


Calculated  and  Observed  Properties  of  Argon 

Calculated  Observed 

X-Alpha  Eicperimental    Experimenlal 

Equilibrium  Equilibrium       Value 

Lattice  Constant            9.66773  10.039          10.039 

Atomic  Mass   t gAr  °       (7.2845784x10^)  (7.2845784x10^)   7.2845784x10^ 

Exchange  Parameter-Alpha     0.7  2131  0.72131            - 

Pressure (in  kilcbars)        -1.21  -3.99            0.0 

Total  Energy              -1C53. 603303  -1053.602791 

Total  Atomic  Energy  ''     -1053.596648  -1053.596648 

Cohesive  Energy             -0,006655  -0,006143       -0.006479 

Zero  Point  Energy          6. 7300x10"  6.2250xl0~'^     5.9600xlO~'^ 

Lennard-Jones  Pararaeters            _A  A               / 

e            7.72920-10  '  7.13456xlO~"    7.55570x10 

a               6.27067  6.51149         6.51149 

Quantum  Parameter  '^       9.03364xl0~  9 .07603x10"  '     8 .57550x10"'^ 
All  quantities  are  in  Rydberg  atomic  units  except  v.liere  noted 
Se-i  references  21  and  36 
See  reference  2j 
See  Appendix  I  for  definition 
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3-2.  The  SCF  Phcncn  Spe.ctta  of  Argon  and  Neon 

As  mentioned  in  section  2-1,  5CH  theory.  Instead  of  conventional  har- 
monic theory,  was  used  to  calculate  the  phonon  spectra  of  argon  and  neon 
because  of  zero  point  motion  considerations.  In  addition,  it  was  felt  that 
SCH  theory  represents  a  conceptual  improvement  over  quasi-harmonic  theory 
for  the  calculation  of  unperturbed  phonons. 

All  phonon  calculations  used  an  extremely  fine  reciprocal  space  mesh 
(large  number  of  phonon  wave  vectors).  Symmetry  considerations  allow  the 
calculation  to  be  done  in  the  irreducible  wedge  of  the  first  BZ  (see 
Figure  A2-1  and  Appendix  II).  Thus,  the  six  non-equivalent  wave  vectors 
used  in  the  selection  rules  and  their  respective  frequencies  represent 
a  small  subset  of  the  calculated  frequencies. 

Due  to  the  three  polarizations  of  the  phonons  and  to  degeneracies  of 
the  phonons  among  these  polarizations,  the  set  of  six  points  corresponds 
to  twelve  energetically  different  phonons.  These  tvrelve  phonons  and  tLeir 
respective  repi-esentatlons  are  listed  as  "Bare  Phonons"  in  Table  3-3  for 
iieon  and  Table  3-4  for  argon. 

The  overbinding  of  neon  in  the  Xa  model  is  reflected  in  the  "bare" 
phoaoT.i  specLra,  This  is  e'.xpected  since,  the  static  solid  cohesive  energy 
is  included  in  the  Lennard-Jones  pair  potential  description.  The  bare 
phonons  for  neon  at  7.7  a.u.  are  a  factor  of  two  higher  than  the  experi- 
m.entally  determined  phonons.  Though  the  8.4348  a..u.  calculation  is  .less 
overbound,  the  "bare"  phonon  fi  equencies  aic  still  substantially  different 
from  the  exuerlmentai  frequencies. 

The  static  lattice  cohesive  energy  (given  in  Tables  3-1  and  3-2)  and 
tlie  average  total  energy  per  particle,  deceii.iLned  in  tlie  d\TiamIc  latt.ice 
using  SCH,  ha>/e  been  u?ed  to  ralcrl  a  to  tlie  7(:'ro   poirit:  energy,  T\\r>    resul- 
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Table  3-3 


Comparison  of  the  Calculated  and  Experimental  Neon  Phonons 


Calculation  Experiment 

a=7.7  a.u,       a=S.434S  a.u.   a=8.4395  a.u. 

Representation       Bare  Perturbed   Bare   Perturbed 
Phonons   Phonons  Phonons  Phonons 


Tis 

0.0" 

0.0 

0,0 

0.0 

0.0 

Ai 

2.327 

2,709 

1.873 

2.194 

1.120 

As 

1.619 

2.063 

1.301 

1.675 

0,805 

xi 

3.375 

3.560 

2.715 

2.873 

1,726 

XI 

2.285 

0      c  -^  *-» 

1.837 

2.045 

1,139 

2i 

2.891 

3.161 

2. 31' 6 

2.559 

1.395 

U 

2.343 

2.544 

1.S84 

2.060 

1.129 

U 

1.512 

2.001 

1.216 

1.618 

0,745 

Ll 

3.389 

3.619 

2.726 

2.927 

1.503 

U 

1.513 

1.974 

1.216 

1.597 

0,752 

wi 

2.267 

2.617 

1.823 

2.110 

- 

W3 

2.886 

3,102 

2.321 

2.505 

- 

Experimental  results  are  from  reference  37  and  are  at  T  =  C.S^K 

SyiTimetry  labels  of  phonons  are  those  C'f  leferencc  29 

c  12 

Frequency  units  are  10    U.z  (THz) 
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Table  3-4 


Compari.soii  of  the  Calculated  and  Experimental  Argon  Phonons 


Calculation  Experiment 

a=9. 66773  a.u.      a=10.039  a.u.   a=10,040  a.u. 


b 

presentation 

Bare 

Perturbed 

Bare 

Perturbed 

Phonons 

Phonons 

Phonons 

Phonons 

Tis 

o.o'^ 

0.0 

0.0 

0.0 

0.0 

Ai 

1.599 

1,848 

1.480 

1.723 

1.342 

As 

1.124 

1.399 

1.040 

1.288 

0.958 

XJ 

2.333 

2.434 

2.159 

2.282 

1.909 

XI 

1.588 

1,762 

1.469 

1.634 

1.363 

Si 

1.994 

2.171 

1.845 

2.021 

1.716 

^3 

1.624 

1.7^5 

1.503 

1.627 

1.356 

U 

1 ,  054 

1.428 

0.976 

1.28A 

0.906 

U 

2.337 

2.486 

2.163 

2.309 

1.915 

LI 

1.055 

1.366 

0.976 

1.270 

0.911 

\Ji 

1.573 

1.814 

1.456 

1.677 

- 

W3 

1.999 

2.143 

1.849 

1.989 

- 

3fi 

Experlniental  results  are  from  reference  38  and  are  for  -.  o^-i^        ^-t  T=10°K. 

Phonons  have  been  scaled  by  (36/AO)"''  " . 

S>'nr-:etry  labels  cf  i>l:onons  sm    from  refere-acc  29 

c  12 

Frequency   \jnits    are    10        i;z    (Tti.'') 
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ting  values  and  the  estimated  trtie  zero   po-fr-t  energy  are  listed  in  Tables 
3-1  and  3-2.  The  overbindirg  cL"  necri  is  se.ii\  agair.  to  be  reflected  in  this 
calci^lated  zero  point  energy. 

The  argon  phonons  cotr.-Dare  more  favorably  with  experiment,  due  to  the 
better  agreement  between  the  calculated  and  experimental  static  cohesive 
energy.  The  calculated  xero  point  energy  for  argon  also  compares  well  with 
experiment. 

3-3.  The  Perturbed  Ph o non  Spe c t r a  of  Argor.  and  Keon 

The  calculation  of  the  perturbed  phonons  is  crude  in  the  sense  that 
a  periodic  macro-cell  of  eight  cubic  unit  cells  (2x2x2)  has  been  used  for 
the  calculiition.  The  A.PW  and  SCK  calculations  were  done  with  larger  peri- 
odic macro-cells  to  ensure  trie  accurate  representation  of  the  solid  based 
on  the  smaller  macro --cell.  This  sma]  J.er  macro-cell  is  equivalent  to  32 
atoc:s  since  the  fee  crystal  structure  has  four  atoms  per  cubic  unit  cell. 
In  the  actual  evaluation  cf  the  EPMc,  expressions  and  selection  rules,  only 
six  non-equivalent  wave  vectors  of  the  electrons  and  phonons  (correspon- 
ding to  32  possible  wave  vectors)  \vere  considered.  As  pointed  cut  in 
Appendix  II,  tliese  six  wave  vectors  have  a.-^sociated  vrith  them  tvrelve  re- 
presentations for  the  plior.'ons  and  the  "p-like"  electrons,  and  six  repre- 
sentations for  the  ''s-lil:e"  electrons. 

Since  tlie  "p-like"'  valence  bands  are  full,  the  transition  from  any  k 
state  in  th-c  p  br.nd  is  only  possible  to  any  empty  k'  state  iii  the  conduc- 
tion bands.  In  view  of  the  ent.rgy  difference  denominator  for  the  second 
order  energy  cc^rrections  (2-42) ,  it  vras  deem.ed  sufficient  to  sum  over  all 
possible  t^'  states  in  the  lov-'e£t  "s-liko"  conduction  band. 

As  the  r,  sum  limit  of  the  /J'W  calculation  v.'as  10,  th.e  E?li]L   expressions 
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of  Appendix  III  could  only  be  suir,med  through  £=9  because  of  the  presence  of 
terras  dependent   on  (5,+l)  .  The  same  square  magnitude  for  the  reciprocal 
lattice  vectors  was  used  as  in  the  APW  calculation. 

The  perturbed  phonons  for  neon  and  argon  are  listed  in  Tables  3-3 
and  3-4  respectively.  Typically  most  bare  phonons  were  perturbed  from 
approximately  ten  to  thirty  percent.  This  is  to  be  expected  in  view  cf  the 
discussion  (see  also  section  2-4)  by  Hedin  and  Lundqvist   of  substantial 
corrections  to  the  bare  phonons  given  by  the  Fr^hlich  Hamiltonian.  However, 
phonons  of  both  neon  and  argon  are  corrected  by  about  the  sane  percentage 
so  the  method  in  general  appears  to  be  reasonable.  That  is,  the  phonon 
corrections  do  not  appear  to  reflect  the  overbinding  situation  in  necn. 

The  neon  perturbed  phonons  are  in  both  cases  poor  because  the  bare 
phonons  reflect  the  overbinding.  The  argon  results  are  substantially 
better,  probably  ss   a  consequence  of  the  good  agreement  between  the  calcu- 
lated and  experimental  static  cohesive  energy.  Though  the  good  agreement 
for  the  argon  cohesive  energy  may  be  fortuitous,  the  method  of  calculating 
the  second  order  energy  corrections  would  seem  to  give  good  results  if  the 
bare  phonons  are  a  reasonably  decent  representation  cf  the  d^TiaTP.ic  lattice 
without  the  electroii-phonon  interaction. 


CBAPTER  IV 


CONCLUSIONS 


The  purpose  of  this  work  has  been  to  develop  a  fessibl.e  first  prin- 
ciples (or  nearly  so)  approach  to  the  calculatiou  of  the  effects  of  the 
electron-phonon  interaction  in  an  insulating  solid.  The  problem  was  forinu- 
lated  in  the  APW-Xa  formalism  as  a  means  of  obtaining  some  computationally 
practical  realization  of  the  calculation.  This  utilization  of  the  Xci  mod&l, 
in  a  way  v;hich  had  not  been  previously  attempted,  clearly  constit:utes  a  te;:.- 
of  the  method.  Nevertheless,  it  was  not  the  purpose  of  this  work  to  w;i,gh 
the  merits  of  various  density  functional  or  local  exchange  schemes.  Rather, 
the  purpose  was  to  see  to  what  extent  the  features  of  a  prototype  density 
functional  model,  with  all  its  admitted  shortcomings,  could  be  exploited 
for  unconventional  purposes. 

The  results  and  discussion  in  section  3--3  indicated  that  the  method 
seemed  to  be  a  reasonable  approach,  even  though  a  rather  crude  calculation 
was  done.  This  is  encouraging  in  view  of  th.e  approximations  made  to 
achieve  the  calcul-Ttion.  Also  it  seem.s  reasonable,  based  on  tlie  results 
of  Chapter  III,  that  better  "uare"  phcnoris  and  crystalline  potentials 
would  give  much  tetter  agreement  v^ith  experiment.  The  unrealisuic  nature 
(af  least  as  comp>-.red  "i  th  the    results  from  potentials  known  to  be  reason- 
able) oi  the  bare  phonous  used  here  can  be  traced,  In  part,  to  the  assumed 
Lennard-'Jones  pair  potential.'^'  At  best,  such  a  pair  potential  is  v.'el3. 
known  to  be  inadequate  to  describe  a  rare  gas  crystal.  Mo;:eover,  we  h^;.ve 
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already  seen  that  the  Xa  incdel  overbinds  peor.  severely  and  argon 
slightly,  so  that  the  calculated  equilibriuni  geometry  is  compressed  vjxth 
respect  to  the  observed  latti:e.   An  ijip roved  fanctional  form  for  the 
electronic  Hamiltonian  would  be  vary  de:3irable  on  these  grounds  alone. 
When  one  considers  the  problems  wlrh  th.e  Xa  band  gap  in  argon,  and  the 
consequent  likelihood  that  tiie  Xa  static  lattice  conduction  band  states 
are  of  less  than  optimum  quality,  the  value  of  an  iaproved  electronic 
Hamiltonian  becomes  even  more  obvious.  In  spits  cf  these  dravzbacks,  it 
remains  a  somewhat  remarkable  fact  that  a  model  as  simple  as  Xci  can  be 
extended  so  far  as  to  yield  a  respectable  estimate  of  phonon  dispersion 
in  a  rare  gas  crystal. 

In  the  present  APV,''  realization  of  the  formailsra,  however,  it  is  not 
practical  computationally  to  undertake  more  refined  calculations.  The 
problem  is  simply  that  of  excessive  compater  time  (see  Appendix  IV).  This 
is  a  direct  consequence  of  the  energy-dependent  APlv  basis,  a  problem  met 
in  otlior  related  work.    A  change  of  basis  functions,  to  say,  a  Gaussian 
type,  (once  the  APW  method  has  been  used  to  calculate  the  eigenvalues  and 
Bloch  functions)  would  reduce  substantially  the  computational  labor  of 
constructing  a  different  basis  set  at  each  eigenvalue.  The  analogue  for 
solids  of  Sambe'r;  molecular  v.^ork,  '   for  instance,  (a  Gaussian  basis  Xa 
method)  would  be  appropriate. 

The  results   represent  an  attempt,  within  certain  approxima- 
tions and  assumptions,  at  a  first  principles  (or  nearly  so)  calculation 
of  microscopic  lattice  dynamics.  A  more  rigorous  approach  would  include 
the  dielectric  function,  and  not  assuir.:-  screening  effects  as  this  work 
has  done.  However,  this  v.'ould  complicate  the  calculation  by  requiring  the 
sddi'"iona]  cor^put  s'^'ion  of  the  dielectiic  function  iir.d    its  inv^'rse,  a 
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difticulc  cask.  Clearly,  the  use  of  thn  Lennard-Jones  form  for  the  unper- 
turbed nuclear  motion  represents  an  intuitively  plausible  circumvention 
of  these  difficulties. 

If  the  method  can  be  refined  and  niade  computationally  practical,  it 
could,  one  might  hope,  be  applied  to  calculations  of  other  effects  caused 
by  the  electroii-phenon  interactioi; .  Based  on  this  work,  the  calculation 
and  methcd  appear,  in  general,  to  be  promising.  This  conclusion  is  proba- 
bly the  most  important  result  cf  this  work:  the  demonstration  of  the  com- 
putational accessibility  of  the  EPI'TE  given  certain  approximations.  With 
the  calculation  now  knovn  to  be  at  least  possible ^  one  may  hope  to  achieve 
better  quantitative  undeirtitandiat^  of  electron-pl'onon  effects  in  a  solid 
by  direct  numerical  investigation. 


APPENDIX  T 
SELF-CONSISTENT  HASI-IONIC  AFPROXIK^TION  aN£  LENNAP.D- JONES  PAIR  POTEInTIALS 

Since  the  unperturbed  or  "bare"  -jhonons  t^ere  calculated  from  cocpu- 
ter  codes  written  by  T.R.  Koehler,  we  shall  follow  the  SCH  treatment  and 
notation  developed  by  Koehler  in  a  series  of  articles.   '  ' '*      Following 
Koehler,  we  direct  our  attention  to  a  one-dimensional  single  particle 
treatment  to  obtain  the  basic  ideas  and  coacepts  of  the  SCH  framework. 
This  treatment  is  only  D.ntended  to  be  represer  tarive. 

Consider  a  particle  v/ith  a  total  Hamiltonian  K  expressed  in  energy 
units  of  c  and  in  length  units  of  C   and  havi-tig  a  m.ass  M, 

H  =  T  +  V  (Al-1) 

2  2  2 

V7here  T  =  -1/2  A  V  and  V  is  the  potential  energy..  A"  is  a  quantum 

7    9  2 

parameter  (=  I'l" /l-io  e) ,    related  to  de  Boer's  parameter'*''  by  a  factor  of  (Ztt)  ' 

Taking  the  hannonic  approximation  to  K  (expanded  about  >~0)  and 

vnriting  it  as  li,  we  get 

h  =  T  +  V°  +  1/2  x^tJ  (Al-2) 

where  V°  ^-  v'   „  and  0  =  V7V  i   ..  The  quantum  mechanical  solution  of 
'x-0  ' x=0 

the  harmonic,  oscillator  is  well  known  ar.d  the  normalized  ground  state 
eigenfunction  is 

|0>  -  (G/T^y^^   exr[-l/2  x^G]  (Al-3) 

,,2    ....  2 
where  &  =  ^/ii  . 
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Taking  the  ground  state  expectation  value  of  R  yields 

E°   =<0iH|0>  =  1/4  A^G  -f  <rj{v!0>  (Al-4) 

where  <Ylwh'>  =  1/2  <'!' [  {VVln^J}  |^''>  has  been  used. 

Minimizing  E°  with  respect  to  G  and  taking  explicit  account  of  the 
normali::ation  gives 

^-''^■^'^-    <olo>   ■     ,oio>2 

To  evaluate  (Al-5),  some  ground  state  matri.;  elerients  are  needed.  Consider 

VlO>  -=■■  -:cGJO>  (Al-6) 

Using  (Al-6)  to  evaluate  'J(i\G>''/v)    yields  the  idenrity 

(  j  0>  V)  1   =  --2<0 1  xGV !  0>  -f  <0  i  v\'  I  0>  (Al-7 ) 

—CO 

However  the  LHS  vanishes  because  the  vrave  function  is  well  behaved  at 
plus  or  minus  infinity;  therefore  we  get 

1/2  G~"^<0!VV!0>  =  <0!xVi0>  (Al-8) 

Operating  on  ((|0>)S)  with  VV  and  using  (Al-7)  and  (Al-8)  gives 

<0|x^'l0>  =1/2  G"^<G|vi0>  -t-  1/4  G"^<0iVVv|0>        (Al-9) 
If  we  let  V  be  unity,  (Al-9)  reduces  to 

<0ix^|0>  -  1/2  g''^<g10>  (Al-10) 


Usi 


sine:  the  identities  (Al-9)  and  (Al-10)  on  (Al-5)  yields 


0  ^  1/4  a'  -  1/4  G"2<^'4i^  (Al-11) 
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which  is  satisfied  provided 


A~G^   =  <0|7WiO>  =  fi-  (Al-12) 


Eqn.  (Al-12)  is  the  self-cons J stent  requirement  and  shows  the 
fundamental  difference  between  ordinary  classical  tiarmonic  theory  and  SGH 
lattice  dynamics.  The  force  constant  is  classical  harmonic  theor}?-  is  the 
second  derivative  of  the  potential  at  equilibrium;  the  force  constant  in 
(AJ.-12)  is  a  ground  state  average  to  be  solved  self-consistently  through 
G.  G  is  seen  to  be  explicit  on  the  LHS  ol  (Al-12)  and  to  be  implicit  en  the 
RHS  through  !o>.  The  extension  of  tlie  above  derivation  to  three  dimensions 
and  N  particles  of  a  crystalline  solid  is  readily  done  (see  Eoehler   ) . 
The  transformation  of  the  problem  to  reciprocal  space  simplifies  the 
actual  numerical  calculation,  just  as  in  ordinary  lattice  d}Tiam.ics. 

The  potential  corresponding  to  (A1-] )  for  the  crystal  problem  is 
assumed  most  commonly  to  be  a  paii-wise  superposition  of  atomic  potentials 

V  =  1/2   l^   v(r..)  (Al~13) 

1,3 

A  fairly      standard   and  v/idely  used  potential    (especially    for   rare 
gas   solids)    is    the  Lennard-Jones    (12-6)    potential 

v(r..)   -   4E[(a/r..)^^   -    (a/r..)^]  (Al-14) 

ij  1.1  13 

T?n's  potential   has   a   zero   at   r..~   o  and   an   energv   tiinimui-i  of   -e  at 
^  13 

r  .  .    -   2  a . 

13 

Since  we  vant  to  fit  e.  cxid   O    to  our  APV7  results,  v.'o  need  the  total 
energy  expression  and  tiu>  volume  derivative  of  ths::  total  energy  expression 
at  equilibrium  volume.  By  choosing  an  N  atom  so3 id  v:ith  our  coordinate 
system  origin  centered  on  an  atom,  vie   may   evaluate  the  total  classical 
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energy  as 

N-1 
E  =  N/2   y   v(r.)  (Al-15) 

1=1    ^ 

where  r   is  the  non-zero  distance  of  the  i   atom  from  the  origin. 
1 

The  stable  lattices  for  argon  and  neon  are  face-centered  cubic  (fee) 
with  the  volume  per  primitive  unit  cell 

U  =  a-^/4  (M-16) 

where   a  is  the  cubic  lattice  parameter.   The  j    atom  in  the  lattice 
is  located  bv 


r  =  a/2  [  m  .e  +  m  .e  +  m.  .e  ]  (M-17) 

-i  xj  X    yj  y    zj  z 

2     2     2   1/2 
with  Ir.i  =  r.  =  a/2  m.  ,   m.  -  (m  ,  +  m  .  +  m  .)     ,  v/ith  m  ,  ,  m   , 

'-3'       .]  3        2        xj       yj       zj  ^-'       yj 

1   -th  , 
m    a  set  of  integers  aporoprxate  to  the  j    atom. 

Using  (Al-IG)  and  (Al-17)  ,  the  following  relations  are  obtained 

r^  =  2  U  m.-^  (Al-lSa) 

J        2 

dr  /dU  -  2/3  m^/r^  (Al-lSb) 

J  2      2 


Miniaiizing   the   total   energy  vrith   respect    to  U 


dE/dU  -    0  -   N/2    2/3   I      nu/r':    dv/aj  .  (Al-19) 

'.2      2  2 

From    (Al-14)    V7e   h.ave 

dv/dr.    -   Ae[-12-^^/r"  +  6//r^]  (Al-20) 

2  2  2 

CombJriing    (Al-16),     (Al-lS)  ,    and    (Al-20)       gives    the    following 

a/a   =    2"^^^    [    I   mT^    /)    if-]'^^  (Al-21) 
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The  lattice  suttis  in  (Al-21)  have  been  evaluated  for  the  fee,  bcc,  sc,  and 
hep  structures.  ■'^''*^''*^  Equation  (Al-21)  shows  that  a/a  is  independent 
of  the  lattice  constant  and  is  determined  solely  by  the  crystal  structure. 
The  cohesive  energy  or  energy  per  particle  may  be  found  using  (Al-15) 
and  (Al-17)  and  solving  for  e 

.  ..  y         1/2  (E/N)    

e  -  Z -ry — r  (Al-2z; 

j  [(a/(m,a/2))-'''  -  (a/im.a/Z)]'] 

These  two  equations,  (Al-22)  and  (Al-21),  together  with  the  appro- 
priate lattice  sums,  have  been  used  to  fit  C   and  a   for  a  representative 
pair  potentir.l  that  best  describes  the  lattice  as  calculated  by  the  M"w 
program.   The  resulting  pair  potential  is  utilized  as  input  by  Koet:ler's 
SCH  progrALi  and  the  "bare"  or  unperturbed  phoiiuns  are  calculated. 


APPENDIX  II 
GROUP  TKEORETICAL  EXPLANATION  OF  THE  TRAN'SITION  SELECTION  RULES 

This  appendix  is  concerned  v.'ith  the  specific  use  and  description  of 
the  group  theory  employed  throughout  this  problem,  particularly  in  deriving 
the  appropriate  selection  rules  and  in  the  labeling  of  the  bands,  phonon 
dispersion  curve;;,  Brillouin  zone,  and  v/ave  vectors.  Also  the  group 
operations  of  all  the  non-equivalent  wave  vectors  and  of  the  uniquf_  vectors 
in  their  corresponding  stars  are  given. 

The  entire  AFW  calculation  was  done  in  the  irreducible  v/edge  of  the 
first  Brillouin  zone  (BZ)  appropriate  to  the  fee  crystal.  The  reduced 
zone  is  shoi-/n  in  Figure  A2-1,  together  with  t'ne  irreducible  wedge  and  the 
two   k_  space  raesh  grids  that  were  used.  The  labeling  convention  of 
Bouci-'aert,  Siroluchowski ,  and  V.'igner  (BSW) '   is  used  throughout  this  disser- 
tation for  the  energy  bands,  phonon  dispersion  curves,  and  all  wave  vector 
labels  (electrons  and  phonons) .  The  list  of  the  six  non-equivalent 
reduced  zone,  vectors  that  were  used  in  deriving  the  selection  rules  is 
given  in  '"able  A2-1, 

The  derivation  of  ihe  selection  rtiJes  is  based  on  the  theory  developed 
by  Zak  '  of  sniall  representation;,  of  the  groups  of  th.e  various  reduced  wave 
vectors.  The  selection  rule  is  calculated  as  the  transition  probability 
between  a  state  of  v;ave  vector   k_  to  a  state  of  wave  vector   k'   b}'  a 
perturbation  or  ei.citation  cf  wave  vectoi  syivjiietry   c|^.   Th.e  actual  calci'- 
lation  is  u.suall.v  done  iu  the  fcllowing  integral  forn 
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Figure  A2-1.  The  first  Brillouin  zone  for  the  fee  structure.  The  points 
labeled  •  are  those  used  in  the  32  point  grid  (which 
corresponds  to  six  non-equivalent  points);  while  those 
points  labeled  O  are  the  256  point  grid  (which  corresponds 
to  nineteen  non-equivalent  points).  The  points  K  and  U  are 
equivalent  and  the  wave  vector  axes  are  in  7r/a  units,  where 
a  is  the  cubic  lattice  constant.  The  irreducible  wedge  (l/48th 
of  the  first  zone)  is  shaded. 
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Table  A2-1 

Representations  and  Non-equivalent  Reduced  Wave  Vectors 

Reduced  Wave  Vector  Representations 

k,  k' ,  ^  ('jiits  of  ir/a)    Phonon/p-like  Valence    s-like  Conduction 

r(o  0  0)                            ri5  Ti 

A(0  1  0)                A,   As  Ai 

X(0  2  0)                Xi  X^  Xi 

Z(l  1  0)               El  Zo  l^  El 

L(l  1  1)                Ll  L§  Li 

W(l  2   0)                Wl  W3  Wl 


'^*,(r)V  (r)f^(r)dr  (A2-1) 

where  ¥,  (r)  is  the  initial  state,  T* ,  (r)  is  the  final  state,  and  V  (r)  is 
k  --  ^.  3.  ~ 

the  perturbation  connecting  the  two  states.  This  integral  vanishes  unless 

k  +  ^=k'  +K(  V7here  K  is  a  reciprocal  lattice  vector)  ,  due  to  lattice 

symmetry. 

Zak's  treatment  11'  done  for  general  non-syniorphic  space  groups  and 

therefore  can  be  somewhat  simplified  for  the  present  work.  We  follow  his 

notation  an.d  convention  and  let  G.  ,  G  ,  G,  ,  stand  for  the  groups  of  the 

k"   q    k 

vectors   k,  q_,  Ic'  respectively,  letting   g  be.  the  number  of  coimnon  group 
elements  and  G  be  the  group  of  coranon  eleiaents  among  the  three  groups 
G  ,  G  ,  Gj^,  ,  one  gets  for  the  selection  rule  formula 

RcG 
c 

where  ^,  .  (R)  ,  E      (?>.),  K-.   i ,  (E-)  are  the  characters  of  the  group  element  R 
ki      q]       k  1 

in  the  i'"  ,  j'"  ,  1"   small  representations  of  k^,  c_,  k'  respectively.  The 
above  formula  actually  only  tells  vliethier  a  certain  transition  is  allowed, 
Hovjever,  by  utilizing  the  vjave  vector  selection  rule  inherent  in  (A2-i) 

k  f  q  -  k'  -f  K  (A2-3) 

we  can  decompose  the  direct  piroduct.s  between  all  sm.all  representations  of 
k  and  £  into  all  allowed  small  reprcsertations  of  k'. 

In  this  work,  we   are  concerned  v.'lth  the  li-ce-centered  cubic  (fee) 
crystal  structure  with  one  at  o;.:  pei  prir.itive  iml  t  cell  (i_._e.  0  (Fm3m) 
space  group  in  the  Schoenflies  and  Into -rrat  ional  Crystallography  notation) 
The  point  group  its-   0   and  in  Table  /.2--2 ,  \;e  give  the  48  group  operations 
of  thib  point  group  fGlloKvim,  the  cviinvention  of  von  der  Lage  and  Eethe. 
Table  A2-3  is  the  group  multiplication  table  lor  the  0   point  group. 
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Table  A2-2 


0,  Point  Group  Ooerations 
h 


Operation  Class 

Number 


0 

^4 


i 


i 


^4 


^4 


^4 


^4 


^4 


10  C, 

4 


11  C2 


Matrix 

Vector 

Representation 

Effect 

10  0 

X 

0  10 

y 

0  0  1 

z 

-1  0  0 

-X 

0-1  0 

~y 

0  0  1 

z 

1  0  0 

X 

0-1  0 

-y 

0  0-1 

-z 

-1  0  0 

-X 

0  10 

y 

0  0-1 

-z 

0-1  0 

-y 

10  0 

X 

0  0  1 

z 

0  10 

y 

-1  0  0 

-X 

0  0  1 

z 

10  0 

X 

0  0-1 

z 

0  10 

y 

10  0 

X 

0  0  1 

z 

0-1  0 

-y 

0  0  1 

z 

0  10 

y 

-10  0 

-X 

0  0-1 

-z 

0  1  0 

y 

1  0  0 

X 

0  ]_  0 

y 

10  0 

X 

0  0- 
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Table  A2-2  coutinued 


Operation  Class  Matrix  Vector 

Number  Representation        Effect 


12  C, 


13  C^ 


14  C2  0-1  0  -y 


15  C2 


16  C2 


17 


18  C^ 


19  C3 


20  C3 


21  C3 


22  C3 


23 


Mat 

rix 

Representation 

0  0 

1 

0-1 

0 

1   C 

0 

-1  0 

0 

0  0 

1 

0   1 

0 

0-1 

0 

-1   0 

0 

0  0 

-1 

0  0- 

-.1. 

0-1 

0 

-1   0 

0 

-1    0 

0 

0   0- 

-1 

0-1 

0 

0  0 

1 

1    0 

0 

0  1 

0 

0   1 

0 

0   0 

1 

1   0 

0 

0   0 

1 

-1    0 

0 

0-1 

0 

0-1 

0 

0   0- 

-1 

1   0 

0 

0   0- 

-1 

-1   0 

0 

0   1 

0 

0-1 

0 

0   0 

1 

-1  c 

0 

0  0- 

1 

1    0 

0 

0-1 

it 

-y 

X 


-y 


-y 


a 
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Table  A2-2  contliiued 


Operation           Class            Matrix  Vector 

Number                       Representation  Effect 

24  ""             C.            0  10  y 

■^            0  0-1  -z 

-10  0  -X 


The  opeiationc  25-43  are  obtained  by  reversing  the  sign  of  the  first  24 
operations;  i.e_.  they  are  the  inversion  operator  (denoted  by  the  class 
symbol   J   )  times  each  of  the  first  2A  operations.  As  an  example,  considei 
operation  4   times  the  inversion:  the  vector  effect  becoines  (x  -y  z), 
the  operation  number  becomes  28,  and  the  class  becomes  JCT . 
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J-!c^^c»^<J•vomcor--oa^r^cNt^<rLnooo^^oa>cMl-^-<rm 

T-i  T-l    iH    tH    fH    iH    ^    ^i    iH    CNI    1— ,    OJ    CN    IM    CN| 


;xi        r-<(Nro<i-mvor~-cxiCr>Oi— (c,  ^n<rm 


fc_»   T— 1   L  •-,   ^  ■  J   ^j    u  t   \^    r^   CO    C^-   O    »— i   CN    C^'J  *cl" 
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r^l  T— I  CNi  r-l 
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o 
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rH 

r^ 
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r-{ 

C^i 
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ao 
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rH 
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VO 

r^ 

CO 
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in 

■■H 

o 

•o 
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a^ 
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CM 
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VO 
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in 
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O 

rH 
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rH 

vD 
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CO 
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in 

CM 

CO 

r^ 

rH 

0^ 

rH 

in 

00 

vD 

o 

rH 

I-- 

r-i 

CO 
CN 

<r 

CM 

o 

C-1 

CNJ 

C-J 
CM 

CO 

CNJ 

CM 

CM 
CM 

CO 
C-J 

CM 

CM 

CO 

o 

CM 

<r 

C^ 

rH 

VD 

CO 

<f 

CM 

O 
C-J 

rH 
CM 

rH 

C-J 

CO 
CM 

-J 
CN 

CO 
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rH 

CM 

-d- 
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.H 
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CM 
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rH 

CM 

rH 

rH 

CO 

i-H 

o 
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rH 
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O-l 
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CN 
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CN 
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CO 

CO 

rH 
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O 

CM 
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o 

CM 

CO 
tH 
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CM 

CM 
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Only  the  filled  p-llke  valence  and  the  unfilled  s -like  conduction 
bands  have  been  considered  and  allowed  to  interact.  The  p-like  representa- 
tions (associated  with  the  vector  k) ,  the  s-like  representations  (associ- 
ated W3.th  the  vector  k'),  and  th=:  representaf: ons  associated  with  the 
phonon  wave  vector  ^  are  listed  in  Tabla  A2-1 .  /J  1  vave  vectors  are  in  the 
reduced  zone  scheme  and  are  in  units  of  T^/a,  where  a  is  the  lattice  coartant. 

In  Table  A2-4,  are  listed  the  group  operations  belonging  to  each 
of  the  32  equivalent  BZ  poinds  used  throughout  chLs  work.  This  table  and 
the  character  tables  of  the  small  represantations  or  tne  wave  vectors 
will  be  used  in  the  illustrative  example  that  will  follow. 

As  an  example  of  the  selection  rule  calculation  consider  an  inter- 
action involving  the  vectors  k  -  (0  1  G)  and  g  -  (1  1  C) ,  and  the  k  repre- 
sei:itation  Ai  and  uhe  q^  representation  T.\.    By  using  (A2-3) ,  all  possible 
allowed  k'  vectors  are  seen  to  be  k'  =  (0  1  0) ,  (1  1  1) ,  and  (1  2  0) . 
Therefore  one  must  use  the  selection  rule  forciulH.  on  all  the  k'  represen- 
tations of  A  (0  1  0),  L  (111),  and  W  (1  2  0). 

Considering  first   only  the  A  representations  of  k ' ,  Table  A2-4  shows 
that  the  groups  of  operations  associated  with  ci  =  (1   1   0)  and  k  =■  (-1  0  C) 

have  only  two  group  elenients  in  common:  the  identity  element  in  class  E 

2 

and  an  iiaproper  rotation  element  in  class  JC,  .   Therefore  the  selectxon 

rule  (A2-2)  reduces  to 

1/2  [S,  (E)e.  (E)C,,(E)  +c       (JC^C.  (JcJ)^,,(JCn]    (A2-4) 
The  character  tables  in  the  KSW  article^^  are  used  to  reduce  (A2-4)  to 

1/2  [^,^(10  +  C.,(dcJ)]  (A2-5) 

which  when  evaluated  over  all  Ai  representations  shows  that  only  the  Ai, 


CO 


Table  A2-4 

Groiip  Operations  of  the  Thirty-two  Equivalent  First  BZ  Points 

Syiumetry   Star  Member  Group  Operations 
Point    (ir/a  units) 

r        (0  0  0)  All  43  Operations 

A        (10  0)  13   7   8  26  28  37  40 

(0  10)  1  4  9  10  26  27  36  39 

(0  0  1)  1  2   5   6  27  28  35  38 

(-10  0)  1  3  7   8  26  2S  37  40 

(0-10)  14  9  10  26  27  36  39 

(0  0-1)  1  2  5  6  27  28  35  38 

X        (2  0  0)  1  2   3  4   7   8  13  16  25  26  27  28  31  32  37  40 

(0  2  0)  1  2  3  4  9  10  12  15  25  26  27  28  33  34  36  39 

(0  0  2)  1  2   3  4  5  6  11  14  25  26  27  28  29  30  35  38 

Z        (1  10)  1  11  2G  38 

(1  -1  0)  1  14  26  35 

(-1  10)  1   14  26  35 

(-1  -1  0)  1  11  26  38 

(1  0  1)  1  12  28  39 

(1  0  -1)  1   15  28  36 

(-1  0  1)  1  15  28  36 

(-]  0  -1)  1      12  28  39 

(0  ]  1)  1   13  27  40 

(0  1  -])  1   16  27  37 

(0  -1  1)  1   16  27  37 

(0  -1  -1)  1   13  27  40 
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Table  A2-4   ccmtinued 

SyriuJierry        Star  Member  Group   Operations 

Point  ("T/a  units) 

I,                   (i   1   1)  1      14   15   16  17  18   25   38   39  40   41   42 

(-1   11)  1     11   12   16  21  22   25   35   36  40  45   46 

(1   -I   1)  1     11   i:   15  19  20  25   35   37  39   43   44 

(-1   -1   1)  1     12   13   14  23  24   25   36   37  38   47   48 

W                   (12   0)                                     1      3  13  16   26   28   31  32 

(2   10)                                     14  12  15   26   27   33  34 

(0  12)                                     14  12  15   26   27   33  34 

(0  2   1)                                     1     2  11  14   27   28   29  30 

(1   0  2)                                     13  13  16   26   28   31  32 

(2   0   1)                                     12  11  14   27   28   29  30 
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A2 5  and  A^  representations  contribute. 

Applying  these  same  techniques  to  the  other  two  cases  (the  L  (1  1  1) 
point  and  VJ  (1  2  0)  point),  gives  the  following  decomposition  between  the 
Aj  and  Ei  representations 

Ai  S  li   =  Ai  6  A2  e  As  ©  Li  Q   Ll  ©  L3  ©  L^  ©  Wi  *  Wl  e  W3   (A2-6) 

where  E  denotes  direct  product  and  <S>   denotes  direct  sum  of  the  represen- 
tations. ^  As  a  check  on  the  decomposition,  the  product  of  the  star  number 
(48  divided  by  the  order  of  the  group)  times  the  dimensionality  of  each 
representation  of  the  LHS  of  (A2-6)  must  be  equal  to  the  sum  of  the  star 
number  times  the  dimensionality  of  each  representation  on  the  RllS .  I'siug 
(A2-6)  as  a.n  exam.ple,  we  see  that  the  LHS  is  equal  to  72  [Aj  dim.  =  1, 
star  =  (48/8)  =  6;  I]  dim,  =  1.  star  =  ('^■S/A;  =  121  and  the  PJIS  is  equal 
to  72  [At  dim.  -  1,  star  ^  6;  A2  dim.  =  1,  star  ^  6;  A5  dim.  ^   2,  star  --^ 
6;  Li  dim.  -  1,  star  =  A;  L2  dim.  -  1,  star  --■  4;  L3  dim.  =  2,  scar  =  4; 
L:'  dim,  =  2,    star  =  4;  Wi  dim.  =  1,  star  =  6;  V,'"2  dim.  =  1,  star  -  6;  W3 
dim.  =  2,  star  -  6].  The  direct  proQuct  is  also  seen  to  be  sj^mnietric,  i_-e_. 
Ai  E  Z,  =  Ii  M  Ai. 

In  tiv?:  actual  calculation  ,  only  the  representations  that  are  in  the 
s--like  excited  ccmdection  band  are  used.  Froir:   (A2-6)   and   Table  A2-I5 
tb.e  only  representations  that  contribute  in  our  example  are  seen  to  be 
A-,  5  Li,  and  I??  •  A  list  of  all  possible  vdence  band  to  conduction  bar:d 
excitations  for  each  q^  (plionoa  vave  vector)  representstJ.on  hati  been  con- 
structed and  is  given  in  Table  A2-5 . 

The  eJectror.-phonon  matrix  element  (EPhK)  calculation  has  been 
simplified  through  the  use  of  group  theory.  Specifically  the  EPillZ  is  first 
calculated  foi"  a  fixed  star  Tiir'n.ber  o*^  t'pe  b '  ^'ector  aad  all  possible  star 


Table  A2-5 
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Table  A2-5   coptinued 
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Table  A2-5  continued 
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Table  A2-5  continued 
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Table  A2-5  continued 


k 

Representation 

Rep: 

3. 

"esentation 

Rep 

/xesentation 

Tis 

B 

Li 

= 

Ll 

Ai 

B 

u 

= 

1:1 

As 

B 

Ll 

= 

2i 

xi 

B 

u 

= 

Ll 

n 

B 

u 

= 

Ll 

xi 

B 

u 

= 

A 

B   Wl 

^3 

B 

u 

- 

Ai 

E, 

B 

u 

=-• 

Nc 

)  Connr.lbution 

Ll 

B 

u 

= 

r 

e  Xi 

hi 

B 

Ll 

= 

Xi 

VI 

B 

Ll 

= 

ll 

W3 

B 

Ll 

= 

l^ 

.  Tis 

B 

u 

= 

Ll 

Ai 

B 

hi 

= 

1:1 

A5 

B 

u 

= 

211 

X4 

B 

Li 

= 

Ll 

x^ 

B 

Ll 

= 

2Li 

2i 

B 

hi 

= 

Ai 

©  v.'l 

J^3 

B 

hi 

= 

A 

e  VI 

X. 

B 

hi 

= 

A 

e  VJ:' 

Ll 

H 

hi 

= 

No 

Cor 

iLribution 

L^ 

B 

hi 

= 

r 

L  e  2x1 

w^ 

B 

LJ 

= 

^1 

V3 

B 

LJ 

= 

^1 

68 


Table  A2--5   continued 
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me^ib^rs  of  the  vector  k,  theii  this  set  of  EPiEs  is  transformed  to  a  Dew 
set  at  esch  of  the  other  members  of  the  star  of  k_'  . 

In  the  computer  program,  the  EPME  is  calculated  for  the  identity 
star  T?_eEber  of  k,'  and  all  the  unique  star  members  of  k.  The  EPME  is  in 
general  a  3x3  matrix  (three  components:  x.,  y,  z;  and  three  polarizations 
or  alternatively  basis  partners)  and  may  be  operated  on  by  the  appropri- 
ate rotation  matrices  of  Table  A2-3.   To  -write  dovm  the  explicit  transfor- 

th 
m.ation,  a  group  operation  of  the  m   star  member  of  the  vector  k  is  denoted 

by  K    and  the  identity  star  member  of  ]c'  and  i    star  member  of  k  of 

the  EP;>r£  is  denoted  by  M-,(k',.,  k,^).  The  transfoi-mation  of  this  EPME  to 

1j.  —  J*  — ^i--- 

the  EPi'iE  of  the  j  '  star  of  ]<_'  and  1    star  of  k  is  accomplished  by 

vith  R  '  the  inverse  operation  to  R  and  the  1  ""  star  member  of  k  deter- 
udncd  by 

R  ,   R,    =  R,  (.A2-8) 

The  necessary  group  cpf")rations  of  the  star  members  of  the  various 
reduced  wave  vectors  are  given  in  Table  A2-5.  This  table  is  based  on  the 
irreducible  v-edge  given  in  Figure  A2-1.  i.e.  all  the  identity  star  members 
of  the  reduced  wave.  \'ectors  are  contained  in  this  wedge. 


BZ 

Star 

Point. 

Generation 

r 

R(0   0   0)    =    (0   0   0) 

A 

R(0  1   0)   =    (0   1   0) 

=    (0  0  1) 

--=    (1   0  0) 

=    (0  -1   0) 
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Table  A2-6 

Star  Group  O7.£rations  Generated  from  the;  S?!x  Mon- Equivalent  BZ  Points 

Group 
Operations 

R  -■  All  -'is  Operations 

R  =  1  4   9  10  26  27   36   39 

R  -  7   13  17  21   32  40  43  47 

R  =  6  11  18  24  29   3S  44  46 

R  =  2   3  12  15   25  28  33  34 

=(0  0  -1)   R  =  8  16   19   23   31  37  41  45 

=  (-1  0  0)   R  =  5   14   20  22  30  35  42  48 
R(0  2  0)=  (0  2  0)    R  =  1  2  3  4  9  10  12  15  2S  26  27  28  33  34  36  39 

=(0  0  2)    R  =  7  8  13  16  17  19  21  24  31  32  37  40  41  43  45  47 

=(2  0  0)    R  =  5  6  11  14  18  20  22  24  25  30  35  33  42  44  46  48 
R(l  1  0)  =  (1  1  0)     R  --=  1   11   26   38 

=  (1  -1  0)    R  =  3  6   28  29 

=  (-1  1  0)    R  =  4   5   27   30 

=  (~.l  -1  0)   R  =  2   14   25   35 

=  (0  1  1)    R  =  10  17   39   43 

=  (0  -1  1)    R  =  12   21   33   47 

=(0  1  -1)    R  =  9   2  3   36   45 

=  (0  -1  -1)   R  -  15  19   34   41 

=  (1  0  1)     r.  -  7   18   40   46 

=  (1  0  -1)   R  -  £  24   37   ^4 

=(-10  3)    R  =  13   2  0   32   48 

=  (-1  0  -x)   R  =  IC   22   31   4  2 
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Table  ^Ji~6   corLtiaued 

BZ  St:  ax  Group 

Point       Generation  Operations 

L     R(l  1  1)  =(111)  R  =  1  14  15  16  17  18  25  33  39  UO   41  42 

=  (-1  1  1)  R  =  3  5  10  13  19  24  27  29  34  37  43  48 

=  (1  -1  i)  R  -  4  6  7  12  22  23  28  30  31  36  46  47 

=  (-1  -1  1)  R  =2  8  9  11  20  21  26  32  33  35  44  45 

W     R(l  2  0)  =  (1  2  0)  R  -  1   3  13  16  26  23  31  32 

=  (2  10)  R  =  5  11  19   21   30  38  41  47 

=(0  12)  R  =  6   14   17   23  29  35  43  45 

=  (0  2  1)  R  =  10  12   22   24   33   39  42   44 

=(10  2)  R  =  2  4  7   S  23   27   3?  40 

=(2  0  1)  R  =  9   15  18  20   34   36  46  48 


APPENDIX  III 

EXx^LIClT  DERIVATION  AND  FORMULAE  FOR  THE  EPME 

We  have  shown  in  Chapter  II,  section  2-5,  that  the  EPME  expression 
to  he  evaluit.ed  uv.nerically  from  (2-63)  using  SAPW  basis  functions  is 

M(S',k';S,K)^:  (^..)'l      II  II      IC*f  ,(K')C«  ^rK).^  ^(R')rf ^(R) 

K'  K  RVX^,  ReG^_   s'  s  ^-  '  "^     "'=-    ^  >^     '^'^ 

X.     X,  XT,     m  ■ —  —  •'       ^ 


'■^  dr    (A3-1) 


£     s'  „   .  Q 


By  considering  an  arbitrary  polarlzaticr  vector  of  the  form 

c   .    =  ae  -f  be  +  ce  (A3-2) 

■■^A     X     y     z 

v:e  can  consider  the  >:,  y,  and  z  components  of  (/O-l)  separately.  This 
is  dor.e  by  rent-mbering  that 

e   -  cosisinOe  -f  sinOsinoe   +  cosBa  (A3-3) 

r  X  y        z 

where  the  angular  convention  of  the  spherical  ham-onics  1^   used.  Taking 

tlie  dot  prorract  between   e    and   C  ,  rives  three  separate  angular  parts, 

r  -'XI'-- 

one    for   each   ccnponent 

Y^     ,(n     .    ,/     ,     .)Y      (0       ,i:      ^i     V*,    ,(e,d)cos4'sinOY„    (e,c>)dO    (A3-4a) 
^,    in        R   k      'K   k'      x.  u.      Rk      Rr/ Ui    !•'    in  i-Ln 


Ws  shall  derive  eiqjlicitly  the  z -component  expression  to  demonstrate 
how  the  final  expressions  were  calculaced.  The  follov>/ing  identities  for  the 
spherical  harTnonics  and  Legendre  pulynoirials  are  based  ou  unpublished  nc^tes 
of  Trickey'   for  use  in  a  k.£  energy  band  method.  The  spherical  hatiaonics 
are  of  the  Condon  and  Shortley  convention   and  are  defined  by 


Y,    (6,6)   =   C„    P„    (cose)exprimd))  (A3-5a) 

Jem  iiti  x^m 

wnere  C.      =    (-J.)      — 7:r^^oT"vrJ  (A3-5b) 

Jem  4.T      \^-^:)  ! 

and  r,   are  the  associated  Legendre  pol\T^.oniiais . 
Jim 

The  identities  that  are  used  in  the  angular  integrals  to  invoke  the 
orthogonality  of  the  spherical  harmonics  are 

cos9fn   =  -7r^^-T^.\r- i„,-,   +  -—- — Y.  ,  J     (A3-6) 


"(21-41)   C,  ,   S",+lm    Cr  ,   £-lm-' 
x-rlm  3c-lm 


.  n^     _       ,.,^,    ^2r.  .(£4t,)(;.+u.-1)v  ( £ -Erf  1 )  ( ^ -lU-f  2 )  „         l//-^7^ 

xra  (2/-rl)    C„  ,   ,     ^-lm-1     C,,  ,-,   ^      »-+lm-l 

sinGY.   --^  exp(-i?)7,,T|f-[— -^ ^^  -L—  y      ]   (A3 -8) 

£-1         ^  (2Mi)  C„,-  , -,  '-^-jm-t-l    C,  T  ,T   £-1ti.!+x 
£+iir.4-A  £-lm+l 

sini  --  1/2 i  (exp[id:]  -  exp[-i<i])      (A3-9) 


cosO  =   1/2  (cxp[ic;;]  +  exp[-i(ii])      (A3-10) 

Y*,.  ,(G,*)Y.  (G,<Odr.'  -  c.,.5   ,  (A3-11) 

1,1 

m=-£ 
where  cosy  --■   u  =  cos9,cos6  +  sinG  sinG  ccs(o  -c,. )  and  1%,  (u)  aie  the 


74 


Legendre  polynomials.    The   greate.iit    siiiir'lifioa*:i.on   of   the   sphericiil   har- 
monic  expressions   is   due    to    the  next    identity  V7l;ich  is   generally  not  well 
known,    but  which   is  provtn   by   Tricka''^      in   his  i-npublished  notes 
dP  dp  £ 

(A3-]3) 

The  LHS  of  (A3-13)  ^&y   be  simplified  by  using  the  following  recursive 
relations  of  the  Legendre  pc  l>Tiouiidls 

dP 

X, 

1-u 


A  =  i^  [  up,,  -  P    ]  (A3-14a) 

du     ,   i      X    x,+l 


d^     -  2      f  ^^f-l  ~  '■''^'J  (A3-14b) 

1-u  ' 

Therefore  (A3-13)  reduces  to 

"^^^-^  -  "  =  ^2-d^  =  r^'    P,(cose.-ucosG^;  +  P^^^vccse^-ucoss^)] 

1-u 

(A3-15) 

It  is  seen  that  ve  have  an  apparent  singularity  et  u  =  ±1  (which  corres- 
ponds to  parallel  or  anti-parallel  directions).,  Hovever,  by  taking  tlie 
appropriate  limit  of  (A3-15)  as  u  -^  ±1,  v/e  may  shox.-  that 

limit  cosO  --—  -   cosQ^— -/"  ""  C±l)  '(£+1)  cos6        (A3-16) 

,,  1   ciu  2   ou  1 

u->---- 1 

By   using    (A3-6)    in    (A3-4c)    and   iirvoking    (A3-11) ,    we   obtain 

^  I     Y*.    ,(e    ,,  ,,6    ,.  ,)Y,    (8   .  ,(*>      )      Y*      ,coseY,    d,Q  - 

x,X      m,n  —  • -2 

;» 

\  J_,^f-(2rS) -(ir^r^         ^mm^^R'h'  -R-k'^^£m^  V 'W 

.  (2.-mrl)  (.f,-hv.-fl)    1/2    .     .„-         ,  ^  ,,         A      ^    %  r^o    -i-7^ 

■*■      ^-CU^-m2£T3)-J         ^^r.^'i^'"K'k'^'£+lm^"Rk'''Rk^    ^  ^^^"^^^ 

The  RH3  of  (A3-17)  may  be  further  simplified  through  the.  use  of 


(A3-13)    .and    (A3-15) 


J(i±ll  r. 


-  ''?,    (cose  ,,  ,-ucose  ,)  +  p,(cose    -ucose  ,,  ,)  )]     (a3-is) 

j?.+i  R'k'  Rk  X  Rk  R'k' 

The  reason  that  the  tei-cis  arp  nor.  coKbined  ir  (A3-18)  is  that  there  are 
Z-dcpandsnt  radial  pares  multiplying  each  terir;  in  zha   parentheses  within 
the  squar?.  brackets.  Coiibining  (A3-18)  and  (A3-1)  gives  for  the  final 
expression  cf  the  z  component 

M  (S',k';n,k)  -   AttI  ^      J        )  [        ">      Jc^"!^    Ct-')c''   (K)r'^   (R')r*''  (R) 

^    "    -        K'  K  R'k  .  ReG  t^    s   ^''-^     -'^  -  n  ,s  -   n,s  - 
—  ^'    ^     k 

[T^      '^Vl^^^^^Rlf'^^^^'^W^  ^  ^£^^°^Vk'-^'^°^V  ' 

■p 

^^Vm^S-^ftN  I   ^Lt  (E'  '^>£^f  (E,r)r^dr 

X-      S   x.-fl     S  ^  O 

For  the  M  a^-^d  K   corponcrts  of  the  EPME,  ve  have  applied  the  above 
X      y    ^ 

methods  and  have  obrained  the  following  e.rpressions 

-    -    -        i'  KR'eG,.,  R^G,^  s-  S   -  '--   n,s-n,s-   n,s- 

I  -7^  f^Vl^^^^^^^Rk-"\k-^-^\'k'«^^^R'k'^  ^ 
£   i  -u  —  —     ^  — 

j,.  fk'R_)j,(kR  )  ^ 


P„(cost^^,^.sinO^^.,,-ucos^:^,^sinGj^.))-^:r^^.^^        u|^^  (E' )fu.  (E)  r^dr 
■     —    —   ii.+±  s'     >'        s   •' o 

■■  ^''m^^^""^R'k'"^"^R'k'-^^"^^"<'Rk"-'^\k^  "■ 
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and 

M    ("',k';n,k)   =   4-;Ti   T      '/        )"  ),'        )      7   f*^^    (K')c'^      (K)r^^      (:R')r*"    (R) 

^                                             K'    K  R'fC,  ,    K-G,     s'    s      '''    '-^  '-'"  ^'    '^ 

—     k     —     k 

V  A^)    [    (p^       (s±n^  ,  3.1r>g      -usin^.    ,,  ,sln5    .,  ,)    + 
'„      ,      2    *■    ^   £+1-          Rk          Fk               rH-'  ?.'k" 

J2,      ].-u  —         —  ■    - 

J         (k'R   )j,(kR   ),^'s  ,  2 

—  ._  j^^  |.]^  S        ■*'  s        c 

-    ^Vl^^'^'R'k^^'^'\'k'--"^'^^*R]^'"W    "" 

^^(-^W^^^k-^^^-^R'k'^-^'k'^^^aF^V':!^!  "i^''^^£+i^"^^'^^-  ^ 

(A3-2i) 


Equafions  (A3 -19),  (A3-20) ,  and  (A3-21)  have  beau  ended  iwro  the 
EPMt  program  (see  Appendix  TA')  togetlie;:  with  the  appropriate  limiting 
expression  (iV3-16)  . 


APPENDIX  IV 
DESCRIPTION  OF  COMPUTER  PROGR^i-IS  AND  CODES 

All  numerical  c-alculations  have  been  done  using  a  set  cf  three 
prograras.  Iwo   of  the  prcgrams  have  been  used  and  tested  in  many  different 
calculations.  The  other  program  was  writtea  hy  the  author  during  the 
course  of  this  calculation. 

One  of  the  established  picgrams,  calcj.lates  the  energy  bands  and 
electronic  wave  functions  based  on  the  API'.'  model,  using  computer  codes 
develope.d  by  Wood  and  othei/^T,  ''    and  vi^ade  self-consistent  and  modified  by 
Connolly  and  others.  ""'''"'" '"  In  the  courbR  of  tb.is  vork,  it  was  found 
necessary  to  mc>dify  the  numerical  precision  (from  single  to  double)  of 
the  card  output  of  the  potential,  eigenvalues,  and  basis  function 
coefficients-  Also  the  basis  function  coefficients  had  to  be  normalized 
properly  foi-  lise  outside  of  the  APU  program  [using  (2-25)]. 

The  oth5:r  established  program  is  Koehler's  SCH  code^   based  on  the 
forr^alisr-i  cutlii^ed  in  Appendix  II.  Tb.is  program  calculates  a  set  of  phonons 
on  any  given  reciprocal  space  mesh  and  at  any  tem,peratu7"e,  given  the 
crystal  st.'uctuT'e,  lattice  para:peter,  lepnard-Joues  parameters,  and  the 
tenipeiature .  This  code  hue   baen  used  and  modified  by  Trickey"'   for  dev;-- 
lopment  of  certain  point  transform  problems. 

The  tiiird  p-rogram.  wass  vritten  entirely  by  tlie  author  except  for 
certain  subroutines  u'liicb  were  borro'.'-ed  from  tbie  /■Ji'V'   program  to  perform 
similiai  tasks,  'fhli^   program  calculates  the  Enir,  over  the  entire  Brillouin 
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zone.  Presently  all  s-like  band  to  p-like  band  transitions  or  all  p-lika 
band  to  s-like  band  transitions  can  be  done,  given  that  one  of  the  bands  is 
full.  The  calculation  of  the  perturbed  phonon  spectrum  is  also  done  \vithin 
this  program. 

It  should  also  be  mentioned  that  the  APW  and  EPMS  progriins  use 
s>Tiiiuetrized  basis  functions.  Group  theory  and  projection  operators  are 
used  to  specify  lists  of  contributing  reciprocal  lattice  vectors  at  any 
given  reduced  zone  vector.  This  "specification  deck"  is  used  as  input  to 
the  APW  and  the  EPMK  programs. 

All  programs  are  written  in  FORTRj:'-N"  anc  have  been  compiled  with 
I.B.K.'s  FOP.TRM  H-Extended  corapilcr  i:al   linVeditcc  and  overlayed  with 
I.B.M.'s  F-128  linkage  editor.   Numerical  precision  is  double  precision 
throughout  all  prograi.is  (approximacely  16  decimal  digits)  . 

The  EPME  program  bar,  been  overlayed  and  arranged  to  run  in  128K  b^/tes 
of  core.  The  progra?;  is  also  designed  to  be  restartable  at  any  point  in  the 
actual  EPHE  calculation.  Th.?.  typical  running  ti.ne  for  the  calculations 
done  in  this  dissertation  was  70  minutes  CPU  (360,/M195)  which  is  approx- 
imately 140  minutes  CPU  (370/M165)  for  all  the  p  to  s  transitions 
on  a  32  point  reciprocal  space  wesh.  The  APVJ  calculations  ran  about  5  to 
10  minutes  CPU  (37Q/M165)   depending  on  the  desired  convergence.  The 
amount  of  required  core  depended  on  the  ba-:is  size.  Ihe  SCH  calculations 
ran  typically  iesr:  than  a  minute  and  in  about  150}.  of  core. 

Num.erical  integration  and  differentiation  are  done  in  the  LPI'LE 
program.  The  potciitia.)  has  to  be  differentiated,  vh^.ch  is  done  by  using 
an  eleven  point  Lagrangian  differentiation  algorithm.  "  The  routine  was 
checked  against  a  J:nov;n  fuiicticn  (  f(r)  =^  l/r  exp[-ar]  )  whose  exact 
derivative  is  calculable.  Th..  numerical  and  exact  derivative  values  agreed 
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to  full  double  precision  for  all  -"-alues  of  r  tested.  This  functior.  was 
chosen  as  a  test  case  because  it  closely  reser.bles  the  numerical  form  of 
the  self-consistent  Ai'W  potential.  The  n-j:?.eri  cal  integration  routine  is 
the  same  as  that  used  in  the  APW  prograTa  and  is  based  on  a  five  point 
Newton-Cotes  integration  algorithn;. '^'*  The  radiaJ  t-7ave  function  of  the 
SchrBdinger  eqn .  is  solved  by  the  Kuiiierov  iriethod-   in  a  subroutine  also 
borrowed  from  the  APVJ  program. 

All  program;  labeling  and  nimbering  con'/entions  of  the  energy/  bands 
and  electron  wave  vectors  follow  BSW"  (see  Appendix  II)  and  are  based  on 
the  "tipecs"  deck  used  for  both  the  /iPU  and  EFME  calculation.  The  same 
notation  is  used  in  J.abeling  the  dispersion  curves  and  phonon  wave  vectors. 
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